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Abstract 

Using 281 ph-^ of e+e collisions recorded at the i/'(3770) resonance with the CLEO-c detector 
at CESR, we determine absolute hadronic branching fractions of charged and neutral D mesons 
using a double tag technique. Among measurements for three and six modes, we obtain 
reference branching fractions B{D^ K-tt+) = (3.891 ± 0.035 ± 0.059 ± 0.035)% and B{D+ 
K~'n~^'K~^) = (9.14 lb 0.10 lb 0.16 ib 0.07)%, where the first uncertainty is statistical, the second 
is all systematic errors other than final state radiation (FSR), and the third is the systematic 
uncertainty due to FSR. We include FSR in these branching fractions by allowing for additional 
unobserved photons in the final state. Using an independent determination of the integrated 
luminosity, we also extract the cross sections (T(e+e~ D^D^) = (3.66 ib 0.03 ± 0.06) nb and 
a{e+e~ D+D-) = (2.91 ib 0.03 ib 0.05) nb at a center of mass energy, E^m = 3774 ib 1 MeV. 

PACS numbers: 13.25.Ft, 14.40.Gx 
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I. INTRODUCTION 



Measurements of absolute hadronic D meson branching fractions play a central role in 
the study of the weak interaction because they serve to normalize many important D meson 
and hence B meson branching fractions. We present absolute measurements of the and 
branching fractions^ for the Cabibbo favored decays K^n^, K^7T~^7r^, 

DO ^ K-7T+71+71-, D+ K-n+7i+, D+ R-tt+tt+ti'^, D+ K'^tt+, D+ K^s^+i:^, 
and — > /T^TT+Tr+Tr", and for the Cabibbo suppressed decay L>+ — > K'^K'n'^. Two 
of these branching fractions, B{D^ — > K'~7r'^) and B{D^ — > ir~7r+7r+), are particularly 
important because most and branching fractions are determined from ratios to one 
of these branching fractions [1]. As a result, almost all branching fractions in the weak 
decay of heavy quarks that involve D'^ oi mesons are ultimately tied to one of these two 
branching fractions, called reference branching fractions in this paper. Furthermore, these 
reference branching fractions are used in many measurements of CKM matrix elements for 
c and h quark decay. 

We previously reported results [2] based on a subset of the data sample used in this 
analysis. The measurements presented here supersede those results. 

We note that the Monte Carlo simulations used in calculating efficiencies in this analysis 
include final state radiation (FSR). Final state radiation reduces yields because D candidates 
can fail the energy selection criteria (the Ai? limits described in Sec. IV) if the energies of 
the FSR photons are large enough. However, many branching fractions used in the Particle 
Data Group (PDG) averages [1] do not take this effect into account. The selection criteria 
imposed in differing analyses correspond to differing maximum photon energies, and hence 
differing FSR effects on the observed yields and branching fractions. Had wc not inchidcd 
FSR in our simulations, our quoted branching fractions would have been lower than we 
report; the difference is mode-dependent, ranging from 0.5% to 3%. 

II. BRANCHING FRACTIONS AND PRODUCTION CROSS SECTIONS 

The data for these measurements were obtained in e^e' collisions at a center-of-mass 
energy Ean = 3.774 GeV, near the peak of the ■(/'(3770) resonance. At this energy, no ad- 
ditional hadrons accompany the D^D^ and D^D~ pairs that are produced. These unique 
DD final states provide a powerful tool for avoiding the most vexing problem in measuring 
absolute D branching fractions at higher energies — the difficulty of accurately determining 
the number of D mesons produced. Following a technique first introduced by the MARK HI 
Collaboration [3, 4], we select "single tag" (ST) events in which either a D or D is recon- 
structed without reference to the other particle, and "double tag" (DT) events in which 
both the D and D are reconstructed. Reconstruction of one particle as a ST serves to tag 
the event as either D^D'^ or D^D~ . Absolute branching fractions for or decays 
can then be obtained from the fraction of ST events that are DT, without needing to know 
independently the integrated luminosity or the total number of DD events produced. 

If CP violation is neghgible, then the branching fractions Bj and Bj for D ^ j and 
D ^ J are equal. However, the efficiencies Cj and for detection of these modes may be 



^ Generally (-D"*") will refer to either D'^ or (D+ or D ), and specification of an explicit D state and 
its decay daughters will imply a corresponding relationship for the D and its daughters. 
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somewhat different since the cross sections for scattering of pions and kaons on the nuclei 
of the detector material depend on the charge of these particles. With the assumption that 
Bj = Bj, the observed yields yi and yj of reconstructed D ^ i and D ^ j ST events will be 

Vi = N^D^iSi and yj = Nj^-^Bjej, (1) 

where A^^d number of DD events (either or D^D~ events) produced in the 

experiment. The DT yield with D ^ i (signal mode) and D ^ j (tagging mode) will be 

Vij = ^DD^i^j^ip (2) 

where ejj is the efficiency for detecting DT events in modes i and j. Hence, the ratio of 
the DT yield (y^j) to the ST yield (yj) provides an absolute measurement of the branching 
fraction Bi, 

B, = (3) 

Due to the high segmentation and large solid angle of the CLEO-c detector and the low 
multiplicities of hadronic D decays, ^ ejCj. Hence, the ratio ej/ejj is insensitive to most 
systematic effects associated with the J decay mode, and a signal branching fraction Bi 
obtained using this procedure is nearly independent of the efficiency of the tagging mode. 
Of course, Bi is sensitive to the signal mode efficiency (e^), whose uncertainties dominate 
the contribution to the systematic error from the efficiencies. 

Finally, the number of DD pairs that were produced is given by 

^nn -f^T7- (4) 
yij cjtj 

Since ~ e^ej, the systematic error for Njj-^ is nearly independent of systematic uncertain- 
ties in the efficiencies. 

Estimating errors and combining measurements using these expressions requires care 
because t/jj and yj are correlated (whether or not i — j) and measurements of Bi using 
different tagging modes j are also correlated. Although D^ and D'^ branching fractions are 
statistically independent, systematic effects introduce significant correlations among them. 
Therefore, we utilize a fitting procedure [5] in which both charged and neutral D meson yields 
are simultaneously fit to determine all of our charged and neutral D branching fractions as 
well as the numbers of charged and neutral DD pairs that were produced (see Sec. IX). The 
input to the branching fraction fit includes both statistical and systematic uncertainties, 
as well as their correlations. We also perform corrections for backgrounds, efficiency, and 
crossfeed among modes directly in the fit, as the sizes of these adjustments depend on 
the fit parameters. Thus, all experimental measurements, such as yields, efficiencies, and 
background branching fractions, are treated in a consistent manner. As indicated above, 
we actually obtain D and D candidate yields separately in order to accommodate possible 
differences in efficiency, but we constrain charge conjugate branching fractions to be equal. 
However, we also search for CP violation by comparing yields for charge conjugate modes 
after subtraction of backgrounds and correction for efficiencies (see Sec. X) . 

We obtain the production cross sections for and D^D~ by combining Nj^q-^q and 

N £)+£)-, which are determined in the branching fraction fit, with a separate measurement of 
the integrated luminosity / Cdt. 



4 



III. THE CLEO-c DETECTOR 



The CLEO-c detector is a modification of tlie CLEO III detector [6-8] in wliicli tlie 
silicon-strip vertex detector lias been replaced with a six-layer vertex drift chamber, whose 
wires are all at small stereo angles to the axis of the chamber [9]. These stereo angles 
allow hit reconstruction in the dimension parallel to the drift chamber axis. The charged 
particle tracking system, consisting of the vertex drift chamber and a 47-layer central drift 
chamber [7], operates in a 1.0 T magnetic field whose direction is along the drift chamber 
axis. The two drift chambers are coaxial, and the electron and positron beams collide at 
small angles to this common axis (see Appendix A). The root-mean-square (rms) momentum 
resolution achieved with the tracking system is approximately 0.6% at p = 1 GeV/c for tracks 
that traverse all layers of the drift chamber. Photons are detected in an electromagnetic 
calorimeter consisting of about 7800 CsI(Tl) crystals [6]. The calorimeter attains an rms 
photon energy resolution of 2.2% at = 1 GeV and 5% at 100 MeV. The solid angle 
coverage for charged and neutral particles in the CLEO-c detector is 93% of Air. 

We utilize two devices to obtain particle identification (PID) information to separate 
from TT^: the central drift chamber, which provides measurements of ionization energy loss 
(dE/dx), and a cylindrical ring-imaging Cherenkov (RICH) detector [8] surrounding the 
central drift chamber. The solid angle of the RICH detector is 80% of Att. As described in 
the next Section, for momenta below 0.7 GeV/c where dE/dx separation is highly efficient 
and RICH separation is not, dE/dx information is used alone. Above this threshold, dE/dx 
and RICH information are combined if both are available. For momenta below 1 GeV/c (the 
entire momentum range of hadrons from D decay at the ijj{3770)) the combined dE/dx and 
RICH particle identification provides excellent separation of kaons and pions, as illustrated 
in Figs. 1 and 2. 
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FIG. 1: Figure (a) shows the efficiency e(vr^) for identifying a pion as a function of the momentum 
p and, on a highly expanded vertical axis scale, (b) shows the probability 7]{K^ — > vr^) that a kaon 
is misidentified as a pion. The reason for behavior of the data above p = 0.7 GeV /c is discussed 
in the text. 



Above 0.7 GeV/c there are modest decreases in e(7r^) and e(i^r^), the efficiencies for 
identifying pions and kaons, respectively; and modest increases in the probabilities, r]{K^ —>■ 
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FIG. 2: Figure (a) shows the efficiency e{K^) for identifying a kaon as a function of the momentum 
p and, on a highly expanded vertical axis scale, (b) shows the probability tj^tt^ — > K^) that a pion 
is misidentified as a kaon. The reason for behavior of the data above p = 0.7 GeV/c is discussed 
in the text. 



TT^) and ri{TT^ K^) of misidentifying a kaon as a pion or vice- versa, respectively. These 
efficiencies and misidentification probabilities are averaged over the whole solid angle of the 
tracking system. However, the RICH solid angle is about 86% of the solid angle of the 
tracking system, and within that solid angle, pion-kaon separation is excellent [8] above 
0.7 GeV/c. Outside of the RICH acceptance, only dE/dx information is available, and 
the lower PID efficiency from dE/dx at high momentum leads to the modest decreases in 
performance observed in this high momentum region. 

The response of the the CLEO-c detector was studied with a detailed GEANT-based [10] 
Monte Carlo (MC) simulation of particle trajectories generated by EvtCen [11], with fi- 
nal state radiation predicted by PHOTOS [12]. Simulated events were reconstructed and 
selected for analysis with the reconstruction programs and selection criteria used for data. 

The integrated luminosity / Cdt (needed only to obtain production cross sections from 
Njyoj)0 and Nd+d-) was measured using the QED processes e"'"e^ e+e~, 77, and fi^fJ.~, 
achieving a relative systematic error of ±1.0%, as described in Appendix C. 



IV. DATA SAMPLE AND EVENT SELECTION 

In this analysis, we utilized a total integrated luminosity of / Cdt = 281 pb^^ of e^e~ 
data collected at center of mass energies near Ecm = 3.774 GeV. The data were produced 
by the Cornell Electron Storage Ring (CESR), a symmetric e^e~ collider, operating in a 
configuration [9] that includes twelve wiggler magnets^ to enhance synchrotron radiation 
damping at energies in the charm threshold region. The rms spread in E^m with the twelve 
wiggler magnets is aE = 2.1 MeV. 



^ The first 56 pb ^ of data were obtained in an earlier configuration of CESR with six wiggler magnets. 
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In each event we reconstructed D and D candidates from combinations of final-state 
particles. Reconstruction begins with standardized requirements for vr^, K^, tt", and Kg 
candidates; these requirements are common to many CLEO-c analyses involving D decays. 

Charged tracks must be well-measured and satisfy track quality criteria, including the 
following requirements: the momentum of the track p must be in the range 50 MeV/c <p< 
2.0 GeV/c; the polar angle 6 must be in the range | cos ^| < 0.93; and at least half of the layers 
traversed by the track must contain a reconstructed hit from that track. Track candidates 
must also be consistent with coming from the interaction region in three dimensions. The 
beams collided close to the origin of the coordinate system, but the collision point in the x-y 
plane (transverse to the axis of the drift chamber system) usually changed somewhat when 
CESR operating conditions changed significantly. Hence, we determined a separate average 
beam position for each data subset bounded by such changes. The period of validity for a 
given average beam position was as short as one run and as long as one hundred runs. (Most 
runs corresponded to a CESR fill and were typically between 40 and 60 minutes long.) For 
each track, we required that the distance d of the track from the average beam position in the 
x-y plane must be less than 0.5 cm {d < 0.5 cm). Finally, we required that the track must 
pass within 5.0 cm of the origin in the z direction {\zo\ < 5.0 cm). The requirements on d 
and zq are approximately five times the standard deviation for the corresponding parameter. 

We identified charged track candidates as pions or kaons using dE/dx and RICH informa- 
tion. In the rare case that no useful information of either sort was available, we utilized the 
track as both a and a tt^ candidate. Otherwise, as described below, we either identified 
it as or tt^, or rejected it if it was inconsistent with both hypotheses. 

If dE/dx information was available, we calculated Xe{'^) Xsi^): where 



from the dE/dx measurements {dE/dx)^^icas, the expected dE/dx {dE / dx)^^.^^ for pions and 
kaons of that momentum, and the measured resolution [a] at that momentum. We rejected 
tracks as kaon candidates when Xe{K) was greater than 9, and similarly for pions. The 
difference x|;(7r) — Xe{K) was also calculated. If dE/dx information was not available, this 
difference was set equal to 0. 

We used RICH information if the track was within the RICH acceptance (| cos^^| < 0.8) 
and its momentum was above 0.7 GeV/c, which is far enough above the Cherenkov threshold 
for kaons that we expect good efficiency for kaons and pions. Furthermore, we required 
that valid RICH information was available for both pion and kaon hypotheses. We then 
rejected tracks as kaon candidates when the number of Cherenkov photons detected for the 
kaon hypothesis was less than three, and similarly for pions. When there were at least three 
photons for each hypothesis, we obtained a difference for the RICH, Xr{'^) ~Xr{K), from 
a hkelihood ratio using the locations of Cherenkov photons and the track parameters [8] . If 
RICH information was not available, we set this difference equal to 0. 

The final particle identification requirement for a kaon (pion) candidate was that the 
track be more consistent with the kaon (pion) hypothesis than the pion (kaon) hypothesis. 
Specifically, we combined the dE/dx and RICH differences in an overall difference, 
Ax^ = x'e{'^) ~ X%{^) + ('?'") ~ Xr(-^)- Kaon candidates were required to have Ax^ > 0, 
and pion candidates were required to have Ax^ < 0. When Ax^ = 0, we utilized the track 
as both a K"^ and a tt^ candidate. 



Xe 
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We formed neutral pion candidates from pairs of photons reconstructed in the calorimeter. 
The showers were required to pass photon quality requirements and to have energies greater 
than 30 MeV. An unconstrained mass M(77) was calculated from the energies and momenta 
of the two photons, under the assumption that the photons originated at the center of the 
detector. This mass was required to be within three standard deviations {3a) of a nominal 
7r° mass value that varied slightly with the total momentum of the tt" candidate. The 
slight change in the nominal 7r° mass compensates for energy leakage in the calorimeter 
for energetic showers. The uncertainty a on Mi^jj) was calculated from the error matrices 
of the two photons; the values of a were typically in the range 5-7 MeV/c^. We then 
performed a kinematic fit of the two photon candidates to the mass M^o from the PDG [1] , 
and the resulting energy and momentum of the 7r° were used for further analysis. 

We built Kg candidates from pairs of intersecting opposite-charge tracks. These tracks 
were not subjected to the track quality or particle identification requirements described 
above. For each pair of tracks, we performed a constrained vertex fit and used the resulting 
track parameters to calculate the invariant mass, M(7r+7r^). We accepted the track pair as 
a Kg candidate if the invariant mass M(7r+7r~) was within 12 MeV/c^ of the mass M^o from 
the PDG [1]. The M(7r+7r") resolution was 2.7 MeV/c^. There is very little background 
under the Kg peak in the M(7r"'"7r~) distribution, so we did not impose requirements on track 
quality or particle identification of the daughters. Also, we did not impose other requirements 
commonly utilized in reconstructing Kg candidates, e.g., requiring that the Kg candidate 
come from the collision point. Imposition of any of these additional requirements would 
have necessitated evaluation of an another systematic uncertainty. 

We formed D and D candidates in the three D° and six D"*" decay modes from combi- 
nations of TT^, K^, 71^ and Kg candidates selected using the requirements described above. 
Two variables reflecting energy and momentum conservation are used to identify valid D 
candidates. 

First, we calculated the energy difference, AE = E — Eq, where E is the total measured 
energy of the particles in the D candidate and Eq is the mean value of the energies of the e"*" 
and beams. The value of Eq was determined from accelerator parameters for each run. 
Candidates were rejected if they failed the IS.E requirements, given in Table 1, which were 
tailored for each individual decay mode. As mentioned in the Introduction, a D candidate 
may be lost if FSR reduces E below the lower limit set by the /S.E requirement. We include 
this effect in our MC simulations. 

Second, we calculated the beam-constrained mass Mbc of the D candidate by substituting 
the beam energy Eq for the energy E of the D candidate, i.e., 

Mlcc' = El-p'c\ (6) 

where p is the measured total momentum of the particles in the D candidate. Valid D 
candidates produce a peak in Mbc at the D mass. To obtain our yields, we fit the Mbc 
distribution for events with Mbc > 1-83 GeV/c^, as described in detail below. 

For the ST analysis, if there was more than one candidate in a particular D or D decay 
mode, we chose the candidate with the smallest \/S.E\. Multiple candidates were very rare 
in some modes, including — > K'n'^ and K^n^n^. and more common in others. 

The largest multiple candidate rate occurred in KgTT~^'K~^'K~ , where approximately 

18% of the events had more than one candidate. 

In two-track events that were consistent with our requirements for D° — > K~7r'^ decays, 
we imposed additional lepton veto requirements to eliminate e'^e" — > e"^e~77, e'^e" — > 
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TABLE I: Requirements on A.E for D candidates. The limits are set at approximately 3 standard 
deviations of the resolution. 



Mode 






Requirement (GeV) 




K 7r+ 




\AE\ < 0.0294 


L) — > 


J\ TT 


TT 


— U.UOoo < l\h/ < U.UooU 




K-TT+ 




|AS| < 0.0200 






-77+ 


< 0.0218 








-0.0518 <AE< 0.0401 


D+ ^ 






|A£^| < 0.0265 


D+ ^ 




-^0 


-0.0455 < A£; < 0.0423 








\AE\ < 0.0265 




K+K 


-7r+ 


A£; < 0.0218 



/i+/i~77, and cosmic ray muon events. We eliminated the event if either the pion or kaon 
candidate track was consistent with being an electron or a muon, utilizing criteria described 
in Appendix B5. A cosmic ray event where the muon has the same momentum as the 
kaon or pion in a decay at rest will peak in Mbc at the beam energy. Removing these 
events simplifies the description of the background shape in the fits. The events from 
e+e~ — > e'^e~77 and e+e" — > ii+ii~j^ populate the Mbc distribution more uniformly. Since 
our DT modes all have at least four charged particles, the e+e~77, fj,+ii~j^, and cosmic ray 
muon event suppression requirements only affect the ST yields. 

In the D+ — > Kg t:+tt+t:~ mode there is a background from Cabibbo suppressed decays 
to D+ — > KgKg7r+. To suppress this background, candidates are rejected if any pair of 
oppositely-charged pions (excluding those from the Kg decay) falls within the range 0.491 < 
M{7r+7r~) < 0.504 GeV/c^. This veto is applied for both ST and DT events. 

To obtain a DT candidate, we applied the appropriate AE requirements from Table I 
to the D candidate and the D candidate in the DT mode. If there was more than one DT 
candidate with a given D and D decay mode, we chose the combination for which the average 
of Mbc(-D) and Mbc(D) — i.e., Mbc = [Mbc(-D) + Mbc(I^)]/2 — was closest to Md- This 
criterion selects the correct combination when an event contains multiple candidates due to 
mispartitioning. (Mispartitioning means that some tracks or 7r°s were assigned to the wrong 
D candidate.) In studies of Monte Carlo events, we demonstrated that this procedure does 
not generate false peaks at the D mass in the Mbc(-D) vs. Mbc(-D) distributions that are 
narrow enough or large enough to be confused with the DT signal. 

V. GENERATION AND STUDY OF MONTE CARLO EVENTS 

We used Monte Carlo simulations to develop the procedures for measuring branching 
fractions and production cross sections, to understand the response of the CLEO-c detector, 
to determine parameters to use in fits for yields, to determine efficiencies for reconstructing 
particular D and D decay modes, and to estimate and understand possible backgrounds. 
In each case e+e~ — > ■^(3770) — > DD events were generated with the EvtCen program [11], 
and the response of the detector to the daughters of the DD decays was simulated with 
GEANT [10]. The EvtCen program includes simulation of initial-state-radiation (ISR) 
events, i.e., events in which the e+ or the e" radiates a photon before the annihilation. 
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The program PHOTOS [12] was used to simulate final state radiation — radiation of pho- 
tons by the charged particles in the final state. We used PHOTOS version 2.15 and enabled 
the option of interference between radiation produced by the various charged particles. FSR 
causes a loss of efficiency due to energy lost to unreconstructed FSR photons; the largest 
effect is a 3% efficiency loss for the decay D° — > X^tt"^. We generated three types of Monte 
Carlo events: 

• generic Monte Carlo events, in which both the D and the D decay with branching fractions 
based on PDG 2004 [13] averages, supplemented with estimates for modes not hsted by 
the PDG, 

• single tag signal Monte Carlo events, in which either the D or the D always decays in 
one of the nine modes measured in this analysis while the D or D, respectively, decays 
generically, and 

• double tag signal Monte Carlo events, in which both the D and the D decay in particular 
modes. 

We applied the same selection criteria for D candidates and DD events when analyzing 
data and Monte Carlo events. We compared many distributions of particle kinematic quan- 
tities in data and Monte Carlo events to assess the accuracy and reliability of the modeling 
of the decay process (event generation) and Monte Carlo simulation of the detector response. 
The agreement between data and Monte Carlo events for both charged and neutral parti- 
cles was excellent for almost all distributions of kinematic variables that we studied. The 
results of this analysis are not sensitive to the modest discrepancies that were observed in a 
few distributions. One exception is the resonant substructure in the multi-body final states 
studied in this analysis. The sensitivity of the analysis to the description of the multi-body 
substructure is discussed further in the section on systematic uncertainties. 

VI. DETERMINATION OF EFFICIENCIES AND DATA YIELDS 

We obtained yields in Monte Carlo events and data with unbinned likelihood fits to 
the distributions of M^c (for single tags) and Mbc(-D) vs. Mbc(-D) (for double tags). We 
determined ST and DT efficiencies from the yields of signal Monte Carlo events. These 
efficiencies include the branching fractions for tt" 'j'j and Kg tt+tt" decays. We 
corrected the MC efficiencies for modes involving Kg daughters to be consistent with the 
updated value of B{Kg — > tt+tt") in the PDG 2006 [1] averages. 

The functions and parameters used to model signals and non-peaking backgrounds in 
these fits are described in Sec. VIA and Appendix A. In Sees. VI B and VIC we discuss 
the fit procedures, the efficiencies, and the data yields for double and single tag events. 
Our procedure was to determine first the parameters describing the momentum resolution 
function in each mode by fitting double tag signal Monte Carlo events where the D and D 
decayed to charge conjugate final states. After determining these parameters, we used them 
in fitting all double and single tag modes in data and Monte Carlo events. 
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A. Signal and Background Shapes and Parameters 



Signal line shapes in the Mbc distributions depend on the beam energy spread, initial 
state radiation from the incident e+ and e~, the ^'(3770) resonance line shape, and momen- 
tum resolution. Appendix A describes the method used to combine these contributions to 
obtain the line shape function that we used to describe signals. 

The Mbc distributions for and events have peaks at Mjjo and M^+j respectively, 
and radiative tails at higher masses due to ISR. The shapes of the peaks are due primarily 
to beam energy spread and momentum resolution. The radiative tails occur at Mbc > Md 
because the momenta of D mesons in events that have lost significant energy due to ISR are 
lower than the momenta of D mesons in events without significant energy loss. Therefore, 
using Eq in Eq. (6) to calculate Mbc leads to Mbc > M^. As described in Appendix A, the 
shape of the radiative tail depends on the resonance line shape and the energy spectrum of 
the ISR photons. 

For the fits to data, our resonance line shape description requires values of the iIj{3770) 
mass and width (M^ and F^, respectively) and the Blatt-Weisskopf radius (r) (see Eqs. (A5) 
and (A6)). The resonance fine shape primarily affects the distribution of the radiative tail 
at Mbc > M^. Hence, our data cannot separate the effects of simultaneous changes to 
the mass, width, and Blatt-Weisskopf radius, and we require external input. The Particle 
Data Group [1] reports three measurements of F^ from MARK I [14], DELCO [15], and 
MARK II [16], of 28 ± 5 MeV, 24 ± 5 MeV, and 24 ± 5 MeV, respectively. The PDG 
averages these to obtain 25.3 ± 2.9 MeV, and it also has a fit which gives 23.0 ± 2.7 MeV. 
Furthermore, there is a recent measurement from EES [17] that gives a width of 28.5 ± 
1.2 ±0.2 MeV. In addition to the width, BES also determines the mass of the ■?/'(3770) to be 
3772.4 ± 0.4 ± 0.3 MeV/c^. In our fits we adopted the BES values for the mass and width^. 
We take the Blatt-Weisskopf radius to be r = 12.3 GeV~^, which is favored by our data 
when M^ and F^ are fixed to the BES values. To assess the systematic uncertainties, we 
vary these parameters as discussed in Section VII. 

We used a sum of three Gaussian functions to describe the momentum resolution of the 
detector. 



G{p; q, dp, fa, Sa, fb, Sb) 



(27r)3/V3 



,)V(2(3a<Tp)2) 



fb 



(SaSbY 



(7) 



Here, q is the true momentum of the D meson; p is its reconstructed momentum; dp is the 
width of the core Gaussian; SaCp is the width of the second Gaussian; /„ is the fraction of 
candidates that are smeared with the width of the second Gaussian; SaSbCp is the width of a 
third Gaussian; and fb is the fraction of candidates that are smeared with the width of the 
third Gaussian. All values of Sa and Sb determined from our fits (see below) are greater than 
2, so the second Gaussian is significantly wider than the first and the third is significantly 
wider than the second. 



^ For the width, we actually used the value 28.6 MeV that appeared in a BES preprint before publication. 
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Combinatorial backgrounds were described by a modified ARGUS function [18] 

a(m;mo,^,p) = Am h - — 2 j ™o/, (8) 

where m is the candidate mass (Mbc), ""^o is the endpoint given by the beam energy, and 
y4 is a normahzation constant. The modification of the original ARGUS function allows the 
power parameter, p, to differ from the nominal value, p = |- The parameters ^ and p were 
determined in each individual ST fit to data or MC simulations. Combinatorial backgrounds 
are very small in DT data, so for DT data and signal MC events, we fixed p = | and used 
values of ^ determined from much larger generic MC samples. 



0140207-006 
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(D°) (GeV/c^) 

FIG. 3: Scatter plot of Mbc(-D) vs. Mbc(-D) for D'^D^ double tag candidates. Signal candidates 
are concentrated at Mbc(-D) = Mbc(-D) = Md. Beam energy smearing ((t(£'o)) smears candidates 
along the Mbc(-D) vs. Mbc{D) diagonal. Initial state radiation (ISR) spreads candidates further 
along the diagonal above the concentration of signal candidates. Detector resolution smears an 
candidate parallel to the Mbc(-D) axis {a{D^)) and parallel to the Mbc{D) axis {a{D^)). Since 
the and resolutions are equal, the resulting distribution is isotropic. Candidates with either 
the or properly reconstructed and the other improperly reconstructed are spread along the 
lines Mbc(-D) = Md or Mbc(-D) = M^. Candidates that are mispartitioned {i.e., where some 
particles are interchanged between the and the D^) are spread along the diagonal. Finally, 
some of the candidates smeared along the diagonal are from continuum events {i.e., annihilations 
to uu, dd, and ss quark pairs) where all particles in the final state are found and used. 

In DT fits, we must include a number of features in our fit function. Figure 3 shows 
the distribution of Mbc(-D) vs. Mbc(-D) for DT D^D'^ event candidates from data, and 
it illustrates the signal and background components in the Mbc(-D)-Mbc(-D) plane. The 
principal features of this two-dimensional distribution are the following. 

• There is an obvious signal peak in the region surrounding Mbc(-D) = Mbc(-D) = Mdo. 
The distribution of the signal candidates in this peak is influenced primarily by beam 
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energy spread, and secondarily by the i/-'(3770) resonance shape and detector resohition. 
The signal also includes a tail due to initial state radiation along the Mbc(-D) vs. Mbc(-D) 
diagonal. This correlation is due to the fact that — neglecting measurement and recon- 
struction errors — the values of Mbc(-D) and M-qc{D) calculated using the beam energy 
will both be too large by the same amount if energy was lost due to ISR. 

• There are horizontal and vertical bands centered at Mbc(-D) = Mdo and Mbc(-D) = M^o, 
respectively. These bands contain DT candidates in which the D (D) candidate was 
reconstructed correctly, but the D (D) was not. 

• There is a diagonal band below the peak that continues through the signal region and 
the radiative tail. This band is populated by the following two sources of background. 

— There are "mispartitioned" DD candidates, in which all of the particles were found 
and reconstructed reasonably accurately, but one or more particles from the D were 
interchanged with corresponding particles from the D {e.g., 7r°s were interchanged 

between the D and the D). 

— There are also continuum events in this band {i.e., annihilations into uu, dd, and 
ss quark pairs). Events fall into this band because all particles in the event were 
reconstructed and used to make the D and the D candidates, so the two candidates 
have equal momentum. 

We accounted for the signal described in the first bullet with the DT signal hne shape 
function given in Eq. (A14). To account for the features in the second and third bullets, we 
included four different background terms in each fit: 

• Two background terms where one of the D mesons is correctly reconstructed and the 
second is incorrectly reconstructed. These terms are described by a signal function of 
Mbc(-D) or Mbc(-D) for the correctly reconstructed D ox D multiplied by an ARGUS 
function of Mbc(-D) or M-qq{D), for the D or D, respectively. 

• One ARGUS background shape in Mbc (defined above) for mispartitioned DD and con- 
tinuum events, multiplied by a Gaussian in AMbc = [-^bc(-D) — Mbc(-D)]/2. The width 
of the Gaussian depends linearly on Mbc- 

• One background term represented by the product of an ARGUS function of Mbc(-D) and 
an ARGUS function of M-qq{D), to account for small combinatorial backgrounds. 

The signal shape parameters describing the effects of detector resolution on the D mass 
(Mbc) distributions are determined by fits to DT signal Monte Carlo samples in which the 
D and D decay to charge conjugate final states. The four parameters controUing the two 
wide Gaussians in the resolution function are then fixed to these values in all other fits, and 
the core resolution dp and D mass values are fixed in all other Monte Carlo fits. The DT 
Monte Carlo samples offer a significantly better signal to background ratio than the single 
tag samples, and there are insufficient statistics to determine these parameters well from 
data. Furthermore, double tag fits allow us to separate the effects of beam energy smearing 
and detector resolution. In single tag fits, the effects of detector resolution and beam energy 
smearing both broaden the Mbc distribution. In double tags, as indicated in Fig. 3, beam 
energy smearing moves the events along the Mbc(-D) = Mbc(-D) diagonal line in a fully 
correlated way while the effects of detector resolution smear events isotropically, including 
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perpendicular to this diagonal. The fitted momentum resolution parameters from Eq. (7) 
are given in Table II. 

TABLE II: The momentum resolution parameters in Eq. (7) obtained from fits to the charge- 
conjugate double tag distributions from signal Monte Carlo events: Gp is the width of the core 
Gaussian, /„ and /& are the fractions of the two wider Gaussians in the resolution function, Gp 
is the width of the second Gaussian, and SaSfe Up is the width of the third Gaussian. 



Mode 






ap (McV/c) 


fa 


fb 


Sa 


Sb 


D^^ 


K-TT+ 




3.73 ±0.13 


0.252 ± 0.040 


0.0081 ± 0.0053 


2.23 ±0.12 


2.92 ± 0.69 


Z?o^ 


K-7r+ 




6.24 ± 0.92 


0.306 ±0.147 


0.0383 ± 0.0146 


2.14 ±0.17 


3.03 ± 0.39 


D^^ 


K-TT+ 


TT~^TT~ 


4.05 ± 0.36 


0.247 ±0.105 


0.0105 ± 0.0050 


2.11 ±0.17 


3.63 ±0.65 


D+ ^ 




-7r+ 


3.95 ±0.22 


0.227 ±0.060 


0.0083 ± 0.0019 


2.16 ±0.10 


4.00 ± 0.24 








4.28 ± 1.44 


0.580 ± 0.170 


0.0498 ± 0.0118 


2.36 ± 0.65 


4.21 ± 0.89 


D+ ^ 




2.13 ±0.75 


0.610 ±0.131 


0.0853 ± 0.0389 


2.49 ± 0.65 


2.24 ±0.21 


D+ ^ 






6.39 ±0.53 


0.300 ± 0.071 


0.0146 ± 0.0132 


2.50 ± 0.24 


3.17 ±2.09 


D+ ^ 






3.69 ±0.62 


0.362 ±0.145 


0.0182 ± 0.0042 


2.16 ±0.19 


5.08 ± 0.85 




K+K 


-7r+ 


4.46 ± 0.21 


0.150 ±0.057 


0.0122 ± 0.0043 


2.12 ±0.20 


3.01 ± 0.45 



B. Double Tag EfRciencies and Data Yields 

We determined double tag yields in data and Monte Carlo events from unbinned maxi- 
mum likehhood fits to Mbc(-D) vs. Mbc(-D) distributions using the signal and background 
functions described in the previous Subsection. The efficiencies, yields from data, and peak- 
ing backgrounds (sec Sec. VII) are given in Tables III and IV for D^D^ and D^D~ events, 
respectively. Since the ARGUS backgrounds are small in signal MC, the errors in the effi- 
ciencies were estimated using binomial statistics. 

The quality of the fits and the small backgrounds in double tag data arc illustrated in 
Figs. 4, 5, and 6. Figure 4 illustrates the Mbc distribution for all DT D^D^ candidates 
combined and for all DT D^D~ candidates combined. The figure emphasizes the fact that 
the DT backgrounds are indeed very small. Figures 5 and 6 illustrate the Mbc distributions 
for DT D candidates in each individual decay mode, tagged with candidates from all of 
the D modes utilized in this analysis. The small peaks visible in the backgrounds in these 
two figures are due to projecting events with a properly reconstructed D and an improperly 
reconstructed D onto the Mbc(-D) axis (see Fig. 3). There are small backgrounds from other 
processes that peak in the signal regions of the Mbc distributions, but are not included in 
the background fit functions. Their sources and how they are handled are described in 
Sec. VII, and their contributions are given in the columns labeled Background in Tables III 
and IV. 

A property of the square-root scales, that are utilized in Figs. 5 and 6, is that all errors 
that are proportional to are the same size on the graph. This results in a better visual 
balance between emphasizing signal (linear scale) or background (logarithmic scale). (The 
error bars for smaller numbers of events are actually somewhat larger than those for larger 
numbers of events because these graphs were plotted with RooFit, and errors in RooFit 
plots are 68% confidence intervals [19].) 
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TABLE III: Double tag efficiencies, yields from data, and peaking background expectations for 
D^D^ events. The efficiencies include the branching fractions for vr*^ — > 77 and Kg vr+vr" 
decays, and the vr*^ and particle identification corrections discussed in Section VIII. The entries in 
the column labeled "Background" are the number of events in the signal peak produced by non- 
signal events and the associated systematic uncertainty; estimation of these values is described in 



Sec. 


VII. The quoted yields 


include these background events. 








Double Ta 


g Mode 




Efficiency (%) 


Data Yield 


Background 






W- 






42.20 ±0.35 


630 ± 25 


< 0.1 




K-TT+ 


W- 


-> K+TT- 




23.11 ±0.31 


1,378 ± 38 


< 0.1 




K-TT+ 


W- 


-> K+TT- 




29.79 ± 0.33 


1,002 ± 32 


11.2 ± 1.6 




K-ir+ir^ 


W- 


-> K+TT- 




23.35 ±0.31 


1,383 ± 38 


< 0.1 






W- 


-> K+TT- 




12.15 ±0.18 


2,679 ± 53 


< 0.1 






W- 


■> K+TT- 




16.17 ±0.27 


1,964 ± 46 


22.1 ± 3.2 






W- 


K+TT- 




30.03 ± 0.33 


955 ± 31 


11.2 ± 1.6 






W- 


K+TT- 


-vrO 


15.97 ±0.27 


1,999 ± 46 


22.1 ± 3.2 






W- 


K+TT- 




20.29 ±0.29 


1,601 ± 41 


33.4 ± 3.4 



3970707-024 



3970707-025 




1.84 



1.86 

/W„.(GeV/c2) 



1.88 



1.84 



1.86 

/Wg_(GeV/c2) 



1.88 



FIG. 4: Projections of double tag candidate masses on the Mbc(-C^) axis for (a) all double tag 
1)0 jjO j^odes and (b) all double tag D~^D~ modes. In each plot, the lines are projections of the 
fit results, the dashed line is the background contribution, and the solid line is the sum of signal 
and background. 



C. Single Tag Efficiencies and Data Yields 

We obtained ST yields in data and Monte Carlo events from simultaneous unbinned 
maximum likelihood fits to the Mbc(-D) and Mbc(-D) distributions for ST D and D events. 
Each fit included a signal line shape function for the signal and an ARGUS function for the 
combinatorial background. 

The signal shape parameters, with the exception of the D mass and the momentum 
resolution ap, were fixed to the values obtained from the fits to the corresponding charge- 
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TABLE IV: Double tag efficiencies, yields from data, and peaking background expectations for 
D+D" events. The efficiencies include the branching fractions for 7r° — 77 and Kg — tt+tt" 
decays, and the vr*^ and particle identification corrections discussed in Section VIII. The entries in 
the column labeled "Background" are the number of events in the signal peak produced by non- 
signal events and the associated systematic uncertainty; estimation of these values is described in 
Sec. VII. The quoted yields include these background events. 







Double Ta 


,g Mode 




Efficiency (%) 


Data Yield 


B ackgr ound 






K~ 7r+ 7T+ 


D~ 


— > K+TT 71 




28.98 


± 


0.33 


2,002 


± 


45 


< 0.1 


D+ 




K~ 7r+ 7T+ 


D~ 


— > K+ 71^ 7l~ 




14.82 


± 


0.26 


685 


± 


27 


< 0.1 


D+ 




K" TT+ TT+' 


D~ 






24.27 




0.30 


272 




17 


4.2 ± 1.1 


D+ 






D~ 


K^Ti-^^ 




13.34 




0.25 


747 




28 


5.8 ± 2.7 


D+ 




K~ 7T+ 7T+ 


D- 




■7r+ 


17.16 


± 


0.27 


404 




20 


8.9 ± 4.3 


D+ 




K~ 7T+ 7T+ 


D' 


K~K+7T' 




24.99 


± 


0.31 


167 


± 


13 


< 0.1 


D+ 




K~7r+7r+7r° 


D~ 


— > K+ 71^ 71^ 




14.90 


± 


0.26 


653 


± 


26 


< 0.1 


D+ 






D~ 


— > K+T^~7^~ 




7.11 




0.20 


213 




17 


< 0.1 


D+ 




K~7X+7X+7X^ 


D" 


s " 




12.18 


± 


0.24 


102 




10 


1.3 ± 0.4 


D+ 




yj-+yj-+yj-0 


D- 






6.15 


± 


0.18 


210 


± 


16 


1.8 ± 0.9 


D+ 




K~ 7T+ 7T+ 7^^ 


D- 




■7r+ 


8.28 


± 


0.20 


125 


± 


12 


2.8 ± 1.3 


D+ 




K"7r+7r+7r° 


D- 


-> K-K+7T' 




12.84 


± 


0.25 


54 


± 


8 


< 0.1 


D+ 




K^gTT+ 


D~ 


— > K+7^~7^~ 




24.29 


± 


0.30 


273 




17 


4.2 ± 1.1 


D+ 




Kl7T+ 


D~ 


— > K+ 7l~ 7^~ 




12.68 




0.24 


102 




10 


1.3 ± 0.4 


D+ 




K^g7l+ 


D- 


s " 




20.55 


± 


0.29 


36 


± 


6 


1.1 ± 0.3 


D+ 




K^s7x+ 


D- 


i^rg 7r-7r° 




10.94 


± 


0.23 


92 


± 


10 


2.1 ± 0.5 


D+ 




Kl7T+ 


D~ 


K'^7T-7r~ 


■7r+ 


15.06 


± 


0.25 


66 


± 


8 


2.0 zb 0.6 


D+ 




Kl7T+ 


D~ 


K~K+7r' 




21.06 




0.29 


23 




5 


0.4 ± 0.1 


D+ 




7r+7r° 


D~ 


— ^ K+7^~7^~ 




13.08 




0.25 


660 




26 


5.8 it 2.7 


D+ 






D~ 


— > K+ 7T~ 




6.29 


± 


0.18 


236 


± 


16 


1.8 ± 0.9 


D+ 






D- 


K9,7T^ 

s 




10.54 


± 


0.22 


94 


± 


10 


2.1 ± 0.5 


D+ 






D- 






5.47 


± 


0.17 


233 


± 


16 


3.8 zb 1.8 


D+ 




7r+7r° 




^ -£V C /( /I 


■7r+ 


7 69 




n 19 

yj.A-ij 


-LOO 


4- 


-LO 


4.0 ± 1.5 


n+ 

U 




jy-O + 


D- 


^ K'K+TT- 




10.61 


± 


0.23 


48 


± 


7 


n r. A_ n n 
yj.o zn u.z 


D+ 




Kl TT+TT+TT- 


D- 


K+Tr~7T~ 




17.30 


± 


0.27 


415 




21 


8.9 zb 4.3 


D+ 




7T+7T+7T~ 


D- 


— > K+7T~7r~ 




8.68 


± 


0.20 


122 




12 


2.8 zb 1.3 


D+ 




Kg 7r"'"7r"'"7r~ 


D- 


K^7T- 




14.75 


± 


0.25 


61 


± 


8 


2.0 zb 0.6 


D+ 






D- 






7.40 


± 


0.19 


136 


± 


12 


4.0 zt 1.5 


D+ 




K'^TT+TT+TT- 


D- 


^K^TT-TT- 


■7r+ 


9.66 




0.21 


87 




10 


3.5 ±1.6 


D+ 




K'^ Tr+TT+TT- 


D- 


K'K+TT 




14.16 




0.25 


33 




6 


0.8 ±0.4 


D+ 




K+K-7r+ 


D- 


K+7T^7T^ 




25.08 




0.31 


169 




13 


< 0.1 


D+ 




K+K-7T+ 


D- 


K+7T~7T~ 




12.47 


± 


0.25 


64 


± 


8 


< 0.1 


D+ 




K+K-7T+ 


D~ 


^K^S^- 




21.17 




0.29 


20 


± 


5 


0.4 ±0.1 


D+ 




K+K-7T+ 


D- 






10.74 


± 


0.23 


76 




9 


0.5 ±0.2 


D+ 




K+K-7r+ 


D- 


^ K^s vr-TT- 


■7r+ 


14.63 


± 


0.25 


39 




7 


0.8 ±0.4 


D+ 




K+K-7T+ 


D- 


^ K-K+TT- 




21.37 


± 


0.29 


13 


± 


4 


< 0.1 



16 




FIG. 5: Projections of double tag D^D^ candidate masses on the Mbc(-C^) axis, with the 
reconstructed in any of the three neutral tag modes. The number of events in each bin is plotted 
on a square-root scale. The lines are projections of the fit results; the dashed line is the back- 
ground contribution and the solid line is the sum of signal and background. Projections of the 
candidate masses on the orthogonal Mbc(-D°) axis are nearly identical to those on the Mbc{D^) 
axis illustrated here. 




FIG. 6: Projections of double tag D^D~ candidate masses on the M-qq{D~^) axis, with the D~ 
reconstructed in any of the six charged tag modes. The number of events in each bin is plotted 
on a square-root scale. The lines are projections of the fit results; the dashed line is the back- 
ground contribution and the solid line is the sum of signal and background. Projections of the 
candidate masses on the orthogonal M-qq{D~) axis are nearly identical to those on the M-qq{D~^) 
axis illustrated here. 



conjugate double tag signal MC samples. The D mass, the momentum resolution cxp, and 
the background ARGUS parameters p and ^ were determined in each of the fits, with the 
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FIG. 7: Distributions of measured Mbc{D) or Mbc(-C) values for single tag and candidates 
with D and D candidates combined in each mode. The points are data and the curves are fits 
to the data. In each plot, the dashed curve shows the contribution of the ARGUS background 
function and the solid curve shows the sum of this background and the signal peak function. The 
number of events in each bin is plotted on a square-root scale. The ST decays are illustrated 
in the left column and the ST Z?+ decays are illustrated in the other two columns. The reference 
modes — > K^ir^ and K~-k~^it'^ are illustrated in the first two plots from the left in the 

top row. 



values of these parameters constrained to be equal for D and D. Figure 7 illustrates the 
Mbc distributions for single tag and data, with D and D distributions combined 
in each plot. Table V gives the ST efficiencies from signal MC and the yields in data. 
These quantities were used in the fit, described in Sec. IX, for branching fractions and 
numbers of DD events. The items in the column labeled Background in Table V are peaking 
backgrounds that were not included in the fit functions; their sources and how they are 
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handled arc described in Sec. VII. 

In the K^n^ mode we note that both the efficiency, determined from Monte Carlo 

simulations, and the data yield are larger for the — > K^7r~ mode than for the — > 
K~7r'^ mode. This is consistent with the larger cross section for hadronic interaction of a 
K" than a K~^. In the CLEO-c detector this manifests itself as a lower particle identification 
efficiency for high momentum K~ than for due to the material in the RICH radiator. 
The difference in the tracking efficiency for K~ vs. is smaller, a few tenths of a percent. 

TABLE V: Single tag efficiencies, yields from data, and peaking background expectations. The 
efficiencies include the branching fractions for 7r° — 77 and K'g tt'^tt" decays, and the 7r° and 
particle identification corrections discussed in Section VIII. The entries in the column labeled 
"Background" are the number of events in the signal peak produced by non-signal events and the 
associated systematic uncertainty; estimation of these values is described in Sec. VII. The quoted 
yields include these background events. 



Sin; 


^le Tag Mode 


Efficiency (%) 


Data Yield 


Background 


DO 


-> K'TT+ 


64.18 ±0.19 


25, 


760 


± 


165 


96 ± 27 




K+TT- 


64.90 ±0.19 


26, 


258 


± 


166 


96 ± 27 






33.46 ± 0.12 


50, 


276 


± 


258 


114 ± 10 




K+TT-TT^ 


33.78 ± 0.12 


50, 


537 


± 


259 


114 ± 10 


DO 




45.27 ±0.16 


39, 


709 


± 


216 


889 ± 135 






45.81 ±0.16 


39, 


606 


± 


216 


889 ± 135 


D+ 


— > K~-K^TT^ 


54.07 ±0.18 


40, 


248 


± 


208 


< 1 


D- 


— ^ if+7r~7r"" 


54.18 ±0.18 


40, 


734 


± 


209 


< 1 


D+ 




26.23 ±0.18 


12, 


844 


± 


153 


< 1 


D- 




26.58 ±0.18 


12, 


756 


± 


153 


< 1 


D+ 




45.98 ±0.18 


5, 


789 


± 


82 


81 ± 22 


D- 




46.07 ±0.18 


5, 


868 


± 


82 


81 ± 22 


D+ 




23.06 ± 0.19 


13, 


275 


± 


157 


113 ± 53 


D- 




22.93 ±0.19 


13, 


126 


± 


155 


113 ± 53 


D+ 




31.70 ±0.24 


8, 


275 


± 


134 


173 ± 83 


D- 


— > Kg TT~TT~TT^ 


31.81 ±0.24 


8, 


285 


± 


134 


173 ± 83 


D+ 




45.86 ± 0.36 


3, 


519 


± 


73 


< 1 


D- 


K-K+TT- 


45.57 ± 0.35 


3, 


501 


± 


73 


< 1 



VII. PEAKING BACKGROUNDS 

In Sec. VIA we described the signal and background shapes used to fit the Mbc distri- 
butions. Monte Carlo simulations indicate that the ARGUS shape used for the background 
provides a good description for the combinatorial background. However, in addition to the 

combinatorial background we also have small backgrounds that peak in the signal region in 
Mbc- These peaking backgrounds arc included in the yields obtained from the fits to Mbc 
distributions, so we must subtract them when we determine branching fractions. In this 
Section we describe what peaking backgrounds we have considered and how we estimate 
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their contributions. Tables III, IV, and V show the background estimates obtained using 
the procedures described below. 

These peaking backgrounds can be categorized as: "internal" backgrounds, where D 
mesons that actually decayed into one signal mode were reconstructed and accepted as 
candidates for a different signal mode; and "external" backgrounds, where decays that we 
do not measure contaminate signals that we are measuring. Our methods for subtracting 
each contribution differ slightly. 

In most cases, we subtracted backgrounds of either type by determining: Bb, the branch- 
ing fraction for a D meson to decay to the background-contributing mode 6; Pb-^i, the 
probability that a D that decays to the mode h is reconstructed as an i candidate; and N^-^, 
the number of D and D mesons produced. We obtain Pb^i from Monte Carlo simulations, 
and Njy^ from the branching fraction fitter. For internal backgrounds we used the values 
of Bb obtained from the fitter, while for external backgrounds, we used fixed values of Bb 
from the PDG 2006 [1] compilation. In practice, at each iteration of the fitter we update 
the background estimates using current Nj^-j^ values, and — for internal backgrounds — the 
current Bb values. This dependence of the subtracted backgrounds on the fit parameters is 
accounted for by the fitter in its minimization. For external backgrounds, we include the 
uncertainties in the PDG values of Bb in our estimates of the systematic errors. 

We identified the major sources of external backgrounds by studying generic DD MC 
samples. We used signal Monte Carlo samples to find the major sources of internal back- 
grounds. In a signal MC event, while the D{D) is forced to decay in a particular signal 
mode, the D{D) decays generically, so external backgrounds will also be present in signal 
MC simulations. To isolate the contributions of internal backgrounds, we removed events in 
which the D[D) decayed in a mode that might contribute external background. 



A. Single Tag Backgrounds 

Doubly Cabibbo suppressed modes (external) Monte Carlo simulations indicate 
that the doubly Cabibbo suppressed decays (DCSD) W K-7r+ and W K-n+n^ 
make the largest contributions to peaking backgrounds for decays to these final states. 
The decay A'^tt+tt+tt^ contributes significantly to the background for that final 

state, but the contribution from the two singly Cabibbo suppressed decays (SCSD), —>■ 
K'K^ Tr+ and K+K^t^-, (see below) is larger. 

However, DCSD should not contribute significant peaking backgrounds to — >■ 
i^~7r"'"7r^ or K^n^n^n^ decays because the charge of the kaon from a DCSD will be 

the same as the total charge of the candidate, so it will not be counted as a signal candidate. 
Hence a particle swap, i.e., a double misidentification — calling a tt"*" a K'^ and the K~ a 
TT" — must also occur for DCSD to contribute peaking backgrounds to these two modes. 
Such particle swaps are quite unlikely and — if they occur — the candidate is less likely to 
satisfy the AE requirement. Monte Carlo studies indicate that DCSD followed by a double 
particle swap should contribute only about one event, so we ignored this background. 

Doubly Cabibbo suppressed decays can contribute to — > Kg tt"^, — > Kg 7r"'"7r'^, and 
— > 7r+7r'^7r~ signals [20], but they are legitimately included in the signals since we 
are measuring decays to Kg rather than K^ or K^ , where Cabibbo favored decays (CFD) 
and DCSD could (in principle) be distinguished. 

The detection efficiency for reconstructing the DCSD — > K^7r~ is the same as it is for 
the signal mode — > K'^n", so Pb-^i — e{K7r) in this case. Hence, the DCSD branching 
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fraction and Nj^q-qo arc all that arc needed for this correction. 

The resonant substructure of — ^ K^tt^h^ is shghtly different for CFD and DCSD 
modes [21], and the same phenomenon is likely to occur for — > i^r+7r~7r+7r~. The dif- 
ferences in resonant substructure may lead to different values of p^^i for these modes. We 
studied this question in MC simulations of these decays, by comparing pb^i for samples gener- 
ated with kinematic distributions fiat in phase space and samples with the nominal Cabibbo 
favored resonant substructure. There is no statistically significant difference between the 
two values of pb^i in either mode. We used the values of pb^i for the fiat distribution when 
estimating backgrounds. 

Z)0 K-K^7z+ and D° K+K^Tr' (external) These SCSD modes can fake 
the decays K^n^n^n^ and K^tt~7T~7i^ , respectively, if the Kg decays to 

71'^ 71'. The probabilities for these backgrounds to appear as signals are suppressed by the 
requirement, described in Sec. IV, that pion tracks originate near the interaction region. We 
did not use an explicit Kg veto to further reduce these contributions. Because we required 
that the pion tracks originated near the interaction region, the Kg momentum spectrum can 
affect Pb^i- We determined Pb^i in signal MC samples for the two decay modes, generated 
with different mixtures of resonant {K*^K^) and non-resonant contributions motivated 
by the PDG averages of previous measurements [1]. There was no statistically significant 
difference in the efficiency for the two mixtures. The factor B{Kg — > tt+tt") is included in 

— multipions (external) Singly Cabibbo suppressed decays can fake decays 
to final states with Kg mesons when a tt'^tt" invariant mass falls within the Kg window. 
We estimated the size of this background by using Kg mass sidebands from data. For the 
sidebands, we required that the reconstructed Kg candidate have a mass in one of the ranges 
0.470 < M{TT+n-) < 0.482 GeV/c^ or 0.5134 < M{ti+ti-) < 0.5254 GeV/c^, and that the 
D'^ candidate using this Kg otherwise satisfied all standard requirements. The Mbc spectra 
of these candidates were then fit with the standard line shapes for the mode being faked. 
The momentum resolutions were set to the values obtained from the charge-conjugate double 
tag fits for these modes in data. 

The yields obtained in the sidebands have a significant contribution from the tails of the 
Kg mass resolution, so some signal is counted in our sidebands. We estimated the magnitude 
of this effect using Monte Carlo simulations and corrected our background estimates. 

Since these background estimates were determined directly from data, they do not depend 
on an input branching fraction or Njy^. 

KgKg 77+ (external) This SCSD mode can be reconstructed as D'^ — > 
KgTi'^TT'^'K'. Wc used two factors to limit this contribution: we vetoed Kg 7r~^7r~^7r~ candi- 
dates in which cither of the tt+tt" combinations satisfied 0.491 < M{7i~^7r~) < 0.504 GeV/c^; 
and we required that the pion tracks originated near the interaction region. 

This final state is dominated by the two-body intermediate state K*~^K^, and thus it is 
modeled well in EvtGen. As the PDG does not fit for this mode's branching fraction, the 
value used for B{D+ K*+K^) was that obtained by the E687 Collaboration [22]. 

Particle swap (internal) A double misidentification — reconstructing a K'^ as a vr"*" 
and a 7r~ as a K~ — can result in a decay being reconstructed as a D° decay. This is 
suppressed relative to correct reconstruction by a factor of ^ 10~^ for D° — > K~7r'^, and is 
not observable in any of the other modes, where the particles have lower momentum and 
better dE/dx discrimination. 

We obtained Pb^i for this process by using the signal Monte Carlo simulations for —>■ 
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K~T[^ . Events with genuine — > K'^t:~ on the other side were rejected, and the yield of 
candidates reconstructed in the remaining events with K^7i~ was measured. 

To verify Monte Carlo simulation of the particle misidentification rate, we reconstructed 
events with two oppositely-charged tracks recoihng against a ^ K'^n", D° — > ir+7r~7r°, 
or — > A'+7r^7r^7r+ tag. These tracks were given particle assignments assuming they con- 
stituted a K^TT^ decay; actual D° — > K~7r~^ events have invariant masses peaking at 
the mass, while D° — > 7t~^71~ and D° — > K~ decays are reconstructed at considerably 
higher and lower masses, respectively. Selecting the events at the mass gives a clean 
sample of decays that are known to be — > K~t^'^, without applying the PID selections. 
We then find what fraction of these events are reconstructed as 7r"'"7r~ and K~K^ after using 
the PID selections, and observe that data and the simulation agree to within 30%. 

Continuum, Radiative Return, and r-pairs We have studied continuum, radiative 
return, and r-pair Monte Carlo samples, and we found no evidence for peaking background 
in any of the signal D decay modes. 

Aside from these backgrounds, there is no indication in the generic DD MC sample for 
other backgrounds exceeding the level of 10~^. 

B. Double Tag Backgrounds 

We calculated double tag background rates separately from single tag rates, by considering 
the same potential sources of background for both the D and D candidates. Because the 
single tag fake rates are small, the probability of a double tag candidate arising from two 
fake single tags was ignored (except as noted below), in comparison to the much higher rate 
from one fake single tag and one real single tag. For the DT background process in which 
D — > i is correctly reconstructed but D ^ k is misreconstructed as a — > J decay, we 
predicted the background event count using 

Ui^k^j = Njj-^ Si Bi Pk^jBk. (9) 

In this equation, p^^j is the probability ioi a D ^ k decay to be reconstructed as a ST 
D ^ J decay. The branching fractions Bi and B^ are taken from the previous CLEO-c 
branching fraction result [2], or the PDG [1] for external modes not included in the earlier 
CLEO-c measurement. Charge conjugate DT backgrounds were set equal. 

An exception to the above procedure occurs for the neutral DCSD modes and the "wrong- 
sign" mode — * K^Kg7r~. Because these fake signals reconstruct as the antiparticle of 
the D that actually generated the signal, it is impossible for them to form part of a double 
tag if the other D was correctly reconstructed. This severely suppresses their contribution 
to double tag backgrounds. We included these decays by choosing a particular wrong-sign 
background mode i, using and Bi as expected for mode i to fake single tags, and then 
summing Eq. (9) over the wrong-sign background modes k for the other side. 

VIII. SYSTEMATIC UNCERTAINTIES 

We take systematic uncertainties into account directly in the branching fraction fit. Ta- 
ble VI and Table VII list the uncertainties that we included in the fit, and a brief description 
of each contribution follows. 
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TABLE VI: Systematic uncertainties and the quantities to which they are apphed in the branch- 
ing fraction fit. Uncertainties not correlated between decay modes are given in the first section, 
and correlated uncertainties in the second. The symbols y and e denote yields and efficiencies, 
respectively. Yield uncertainties are additive and efficiency uncertainties are multiplicative. See 
the text for the distinction between e(Charged) and e{K^). The detector simulation uncertainties 
are determined per charged track or per neutral pion or kaon. Uncertainties for other efficiencies 
are determined per D. In addition to the systematic uncertainties listed here, we apply five more 
mode-dependent systematic uncertainties listed in Table VII. 



Source 


Uncertainty (%) 


Quantity or Decay Mode 


DT Signal Shape 


0.2 


y(An DT Modes) 


Double DCSD Interference 


0.8 


y (Neutral DT) 


Detector Simulation 


0.3 


e(Charged) Tracking 




0.6 


e{K^) Tracking 




1.8 






2.0 


6(^0) 




0.25 


e(7r±) PID 




0.3 


e{K^) PID 


Lepton Veto 


0.1 


e{D^ K-TT+) ST 


Trigger Simulation 


0.2 






0.1 


e(L>+ ->i^0 7r+) 


|A£^| Requirement 


1.0 


e{D+ K'^TT+TT^) and e(L>+ ^ K+K-n+) 




0.5 


e(An Other Modes) 



TABLE VII: Mode-dependent systematic uncertainties. The systematic uncertainties for the signal 
shapes are correlated among all ST modes. The systematic uncertainties for FSR are correlated 
among all ST and DT modes. Other uncertainties are uncorrelated. The background and sig- 
nal shape uncertainties are uncertainties on the yields, the other uncertainties in the table are 
uncertainties on the efficiency. Yield uncertainties are additive and efficiency uncertainties are 
multiplicative. 



Mode 


Background 
Shape (%) 


ST Signal 
Shape (%) 


FSR (%) 


Resonant 
Substructure (%) 


Multiple 
Candidates (%) 




0.4 


0.3 


0.9 




0.0 




1.0 


0.5 


0.3 


0.3 


0.8 


D° ^ K-TT+TT+TT- 


0.4 


0.7 


0.8 


1.2 


0.0 


D+ K-'K+'K+ 


0.4 


0.3 


0.7 


0.6 


0.0 


D+ K-TT+TT+TT^ 


1.5 


1.3 


0.3 


0.5 


0.5 


D+ 7r+ 


0.4 


0.4 


0.5 




0.2 




1.0 


0.5 


0.1 


1.2 


0.0 


D+ K^TT+TT+TT- 


1.0 


0.6 


0.6 


0.5 


0.0 


D+ K+K-7T+ 


1.0 


0.6 


0.3 


1.3 


0.2 



Signal shape (DT and ST) We gauge the sensitivity of the ST and DT yields to 
variations in the Mbc fit functions by repeating the fits with alternative fit functions. We 
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vary the parameter values of the signal line shape. The main parameters here are the width 
and mass of the ip{3770) as well as the Blatt-Weisskopf radius. We vary these parameters 
by ±2.5 MeV, ±0.5 MeV/c^, and ±4 GeV~^ respectively and combine the changes in the 
yields in quadrature to obtain the systematic uncertainty assigned in each mode. We also 
vary the resolution function parameters /„, fb, Sa, and Si,. 

We have also tried alternative forms for the parametrization of the line shape; in particular 
we tried using the form used by Mark II [23]. The event yields we determine are insensitive 
to the parameterization of the line shape. 

Double DCSD interference In the neutral DT modes, the CFD amplitudes can inter- 
fere with amplitudes where both and undergo DCSD. This interference is controlled 
by the DCSD/CFD rate ratios (Rws) and relative phases {6). If we assume common values 
of Rws and 6 for the three modes, then the relative size of the interference effect is 
A 2Rwscos25. Because of uncertainties in the value of 5, we assign yield uncertainties 
of 0.8% to span the allowed range of A for Rws — 0.004, which is approximately the mea- 
sured value for D — Kit [1]. These conservative uncertainties are applied incoherently to 
all neutral DT yields. 

Detector simulation — Tracking and Kg efficiencies We estimate uncertainties due 
to differences between efficiencies in data and those estimated in Monte Carlo simulations 

using the partial reconstruction technique described in Appendix B. No significant biases 
are found. A tracking efficiency systematic uncertainty e(Charged) of 0.3% is applied to each 
K"^ candidate and each tt"^ candidate (including those from Kg — > tt+tt" decay). This uncer- 
tainty is fully correlated among all charged tracks in the event. An additional 0.6% tracking 
systematic uncertainty e{K^) is applied to each K"^ track; this uncertainty is not correlated 
with the 0.3% uncertainty for all charged tracks, but it is correlated among all charged kaons. 
The charged kaon systematic contribution arises from a two-standard-deviation discrepancy 
between data and MC simulations in the relative K'^ and K" efficiencies (Appendix B). To 
be conservative, we have assigned this additional uncertainty even though we find no such 
discrepancy in the relative 7r+ and vr" efficiencies, or in the average and the average n'^ 
efficiencies. A Kg reconstruction efficiency systematic uncertainty of 1.8% is applied to Kg 
candidates, correlated among Kg candidates. 

Detector simulation — tt" efficiency Possible differences in 7r° reconstruction effi- 
ciency between data and Monte Carlo simulations are also investigated using the partial 
reconstruction technique described in Appendix B. We find a small bias and correct for it 
by multiplying the efficiencies determined in Monte Carlo simulations by 0.961", where n 
is the number of reconstructed tt^s in each final state. The efficiencies fisted in Tables III, 
IV, and V include this correction. We assign a correlated systematic uncertainty of 2.0% to 
each 7r°. 

Detector simulation — Particle Identification efficiencies Particle identification 
efficiencies are studied by reconstructing decays with unambiguous particle content, such 
as D° X^TT+TT- and K+K'. We also use K'tt+t:^, where the K' and 

TT"*" are distinguished kinematically. The efficiencies in data are well-modeled by the Monte 
Carlo simulation with small biases. We correct for these biases by multiplying the efficiencies 
determined in Monte Carlo simulations by 0.995' x 0.990™, where / and m are the numbers of 
PID-identified tt^s and PID-identified K^s, respectively, in each final state. The efficiencies 
listed in Tables III, IV, and V include these corrections. We assign correlated uncertainties 
of 0.25% and 0.3% to each tt^ and K^^ respectively. We do not assign these corrections and 
uncertainties to K^ daughters, because they are not subjected to the tt^ PID requirements. 
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Lepton veto As discussed in Section IV, in events with only two tracks we required 
— » K^n^ ST candidates to pass additional requirements to eliminate e+e~ e"'"e~77, 
e~^e~ — > jj,'^ iJ,~ 'f'y , and cosmic ray muon events. These requirements eliminate approximately 
0.1% of the real — > K'n'^ candidates, and we include a systematic uncertainty of 0.1% 
to K^n'^ ST yields to account for the effect of these additional requirements. 

Trigger simulation Most modes are efficiently triggered by a two-track trigger. How- 
ever, in the modes — K~Tr~^7:^ and —>■ Kg7r~^, Monte Carlo simulation predicts 
a small inefficiency (0.1-0.2%) because the track momenta may be too low to satisfy the 
trigger or because the Kg daughter tracks may be too far displaced from the interaction 
region. For these two modes, we assign a relative uncertainty in the detection efficiency of 
the size of the trigger inefficiency predicted by the simulation. 

|A£?| requirement Discrepancies in detector resolution between data and Monte Carlo 
simulations can produce differences in the efficiencies of the AE requirement between data 
and Monte Carlo events. No evidence for such discrepancies has been found, and we include 
systematic uncertainties of 1.0% for K^^tt^h^ and K^K~7r^ decays, and 0.5% 

for all other modes. These uncertainties are taken to be correlated in e^, e^j, and e^j. 

Background shape We estimate the uncertainty in ST yields due to the background 
shape by repeating the ST fits with alternative background shape parameters. These al- 
ternative parameters are determined from the Mbc distributions of events in high and low 
AE sidebands. For each mode, we fit each sideband with an ARGUS function to determine 
shape parameters and then repeat the ST yield fits with the ARGUS parameters fixed to 
these values. The resulting shifts in the ST yields are used to set the value of the systematic 
for each mode. 

Final state radiation In Monte Carlo simulations, the reduction of DT efficiencies due 
to FSR is approximately a factor of two larger than the reduction of ST efficiencies due 
to FSR. This leads to branching fraction values larger by 0.5% to 3% than they would be 
without including FSR in the Monte Carlo simulations. We assign conservative uncertainties 
of ±30% of the FSR correction to the efficiency as the uncertainty in each mode. This 
uncertainty is correlated across all modes. 

Resonant substructure The observed resonant substructures of three- and four-body 
decay modes in our simulations are found not to provide a perfect description of the data. 
Such disagreements can lead to wrong estimates of the efficiency in the simulation. We 
estimate systematic uncertainties for the three- and four-body modes from the observed 
discrepancies. These uncertainties in efficiency are not correlated between modes, but the 
correlations in systematic uncertainties for the efficiency of mode i are taken into account 
in ei, eij, and eu. 

Multiple candidates In our event selection, we chose a single candidate per event per 
mode. So, in general, because the correct candidate was not always chosen, our signal 
efficiencies depend on the rate at which events with multiple candidates occur. Using signal 
Monte Carlo samples, we estimate the probability of choosing the wrong candidate, V, when 
there are multiple candidates present. We also study the accuracy with which the Monte 
Carlo simulations model the multiple candidate rate, TZ, in data. If V is non-zero and if TZ 
differs between data and Monte Carlo events, then the signal efficiencies measured in Monte 
Carlo simulations are systematically biased; if only one of these conditions is true, then 
there is no efficiency bias. Based on the measured values of 7'(7^data/'^MC ~ 1); we assign 
the systematic uncertainties shown in Table VII to ST efficiencies. For each decay mode the 
multiple candidate systematic is correlated between the D and D decay for single tags. 
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Luminosity For the e+e~ DD peak cross section measurements, we include additional 
uncertainties from the luminosity measurement (1.0%). The luminosity measurement and 
the uncertainties are discussed in detail in Appendix C. 

IX. BRANCHING FRACTION FITS 

To determine the nine branching fractions as well as Njjo-^o and Nd+o-, we perform 
a single fit that takes as input our measured event yields and efficiencies for the 9 ST 
modes and 45 DT modes given in Tables 111, IV, and V. In this branching fraction fit, we 
correct these event yields not only for efficiency but also for crossfeed among the ST and DT 
modes and for backgrounds from other D decays. The estimated crossfeed and background 
contributions induce yield adjustments of no more than 4%. The dependence of these 
adjustments on the fit parameters is taken into account both in the yield subtraction and in 
the minimization. In addition to the correlated and uncorrelated systematic uncertainties, 
the statistical uncertainties on the yields, efficiencies, and background branching fractions 
are also included in the fit. 

We validated the algorithm and the performance of the branching fraction fit — as well 
as our entire analysis procedure — by measuring the branching fractions in generic Monte 
Carlo events. We find that the results of this procedure are in excellent agreement with the 
input branching fractions used in generating the events; the measured branching fractions 
and DD yields were all within 1.5 standard deviations of the input values. The overall 

of the difference between the fit results and the Monte Carlo inputs, accounting for the 
correlations among the fit parameters, is 13.6 for 11 degrees of freedom, corresponding to 
a confidence level of 26%. Furthermore, the generic Monte Carlo sample has an order of 
magnitude more events than our data, so the statistical errors in this test are about a factor 
of three smaller than in data. The systematic uncertainties are also substantially smaller 
than those estimated for data, so the agreement between measured and generated branching 
fractions of the generic Monte Carlo events is a stringent test of our entire analysis procedure. 

The results of the fit to data are shown in Table VIIL The of the fit is 39.2 for 52 
degrees of freedom, corresponding to a confidence level of 98%. The fit, which includes 
statistical and systematic errors for the input measurements yields total errors for the fit 
parameters. The statistical errors for these parameters are determined separately from a fit 
that includes only the statistical errors of the inputs to the fit. Then the systematic errors 
are determined from the quadrature differences between the total errors and the statistical 
errors. We also repeat the fit after removing the FSR systematic uncertainties for the 
efficiencies to obtain the separate contributions of the FSR uncertainties to the systematic 
errors. If no FSR had been included in the simulations to calculate signal efficiencies, then 
all of the branching fractions would be 0.5% to 3% lower. We hst the shift ApsR for each 
mode in Table VIII . 

Table IX gives the correlation matrix for the eleven fit parameters. In the absence of 
systematic uncertainties, there would be no correlation between the charged and neutral D 
parameters. 

The ratios of branching fractions to the reference branching fractions given in Table VIII 
are not free parameters in the fit, but are derived from the fitted branching fractions. These 
branching ratios have higher precision than the constituent branching fractions. The total 
errors (statistical and systematic) are calculated using the correlation matrix in Table IX. 
Statistical errors for the branching ratios are obtained using the correlation matrix derived 
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TABLE VIII: Fitted branching fractions and DD pair yields. For Nj^q-^q and Nd+j^,^, uncertain- 
ties are statistical and systematic, respectively. For branching fractions and ratios, the systematic 
uncertainties are divided into the contribution from FSR (third uncertainty) and all others com- 
bined (second uncertainty). The column of fractional systematic errors combines all systematic 
errors, including FSR. The last column, ApsR, is the relative shift in the fit results when FSR is 
not included in the Monte Carlo simulations used to determine efficiencies. 



Parameter 



Fitted Value Fractional Error ApsR 

Slal.(S{) Sysl.(Si:) (%) ' 



0.8 


1.3 


+0.1 


0.9 


1.8 


-3.0 


0.8 


2.7 


-1.1 


0.9 


2.4 


-2.4 


1.0 


1.2 


+0.1 


1.1 


1.9 


-2.3 


1.3 


2.8 


-1.0 


1.4 


2.5 


-1.8 


1.3 


3.5 


-0.4 


1.5 


3.0 


-1.9 


1.8 


2.6 


-1.2 



B{D^ K-tt+ttO) 

b\d^ K-TT+TT+TT-) 

Nd+d- 

B{D+ K-tt+tt+ttO) 
B{D+^K^7T+) 

B{D+ ^ K^TT+TT^) 

B{D+ Kg7r+7r+7r-) 

B{D+ -> K+K-TT+) 



(1.031 + 0.008 + 0.013) X 10^ 
(3.891 + 0.035 + 0.059 + 0.035)% 
(14.57 + 0.12 + 0.38 + 0.05)% 
(8.30 + 0.07 + 0.19 + 0.07)% 
(0.819 + 0.008 + 0.010) X 10^ 
(9.14 + 0.10 + 0.16 + 0.07)% 
(5.98 + 0.08 + 0.16 + 0.02)% 
(1.526 + 0.022 + 0.037 + 0.009)% 

(6.99 + 0.09 + 0.25 + 0.01)% 
(3.122 + 0.046 + 0.094 + 0.019)% 
(0.935 + 0.017 + 0.024 + 0.003)% 



B{D^ K-Tr+TT^)/B{K-TT+) 

B{D^ K-TT+TT+TT-)/B{K-TT+) 

B{D+ K-7r+7r+7rO)/'S(K-7r+7r+ 

b\d+ K0 7r+)/B(K-7r+7r+) 

B{D+ K0 7r+7rO)/i3(K-7r+7r+) 

B{D+ Kg7r+7r+7r-)/5(K-7r+7r^ 

B{D+ K+K-Tr+)/B(K-Tr+Tr+) 



3.744 + 0.022 + 0.093 + 0.021 
2.133 + 0.013 + 0.037 + 0.002 
0.654 + 0.006 + 0.018 + 0.003 
0.1668 + 0.0018 + 0.0038 + 0.0003 

0.764 + 0.007 + 0.027 + 0.005 
0.3414 + 0.0039 + 0.0093 + 0.0004 
0.1022 + 0.0015 + 0.0022 + 0.0004 



0.6 2.6 +1.9 

0.6 1.7 +0.5 

0.9 2.7 +1.4 

1.1 2.3 +0.5 

0.9 3.5 +2.0 

1.1 2.7 +0.4 

1.5 2.2 +1.1 



with only statistical errors. The systematic errors are then obtained from the quadrature 
difference between the total and statistical errors. These branching ratios are also sensitive 
to final state radiation, and — without these corrections — all would be 0.5% to 2% higher. 

We obtain the e+e~ — > DD cross sections by dividing the fitted values of Nj^q-^q and 
Nd+d- by the collected luminosity, / Cdt = 281.5 ± 2.8 pb""*^ (see Appendix C). Thus, at 
Ecm = 3774 ± 1 MeV, we find the values of the production cross sections given in Table X. 
(The uncertainty of 1 MeV corresponds to the range of center-of-mass energies in our data 
sample.) 



X. CP ASYMMETRIES 



Although this analysis assumes equal rates for decays to charge-conjugate final states / 
and /, the separately determined yields and efficiencies for charge-conjugate decays allow 
us to calculate CP asymmetries, 

AcAI) - (10) 
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TABLE IX: The correlation matrix, including systematic uncertainties, for the branching fractions 
and numbers of DD events determined from the fit. 





N^Q-^o Ktt Ktttt^ K'k-k'k 


Nd+d- Kit-k Kit-kit^ K^-k K^-kit^ KItttt-k KKtt 


N — 

B{K-Ti+) 

^(i^-7r+7r°) 

B{K- 7^+1^-7^+) 


1 -0.65 -0.34 -0.41 
1 0.44 0.70 
1 0.38 
1 


0.39 -0.19 0.01 -0.14 -0.09 -0.08 -0.09 
-0.22 0.52 0.23 0.28 0.15 0.30 0.35 
-0.11 0.28 0.66 0.14 0.51 0.17 0.21 
-0.09 0.51 0.29 0.28 0.17 0.37 0.34 


Nd+d- 

B{K-7T+7T+) 
B{K-TT+TT+TT^) 

BiKl7r+) 

^(/4^+7r+7r-) 

B{K+K-7v+) 




1 -0.61 -0.24 -0.48 -0.30 -0.33 -0.38 
1 0.43 0.52 0.32 0.51 0.55 
1 0.27 0.56 0.29 0.32 
1 0.55 0.72 0.31 
1 0.50 0.20 
1 0.30 
1 



TABLE X: Production cross sections for e'^e^ DD and the ratio of D^D^ to D^D^ cross 
sections. The uncertainties are statistical and systematic, respectively. The charged and neutral 
cross sections have a correlation coefficient of 0.57 stemming from systematic uncertainties and 
from the common use of the luminosity measurement. 



Quantity Value 

a{e+e- D^W) (3.66 ± 0.03 ± 0.06) nb 

C7(e+e- D+D-) (2.91 ± 0.03 ± 0.05) nb 

(7(e+e- ^ DD) (6.57 ± 0.04 ± 0.10) nb 



(7(e+e- ^ £)+D-) 
a(e+e- ^ D^W) 



0.79 ±0.01 ±0.01 



for each mode /. In this expression, the CP asymmetry Acp{f) is calculated from n{f) and 
n(/), the single tag yields obtained for the charge conjugate modes / and /, after subtraction 
of backgrounds and correction for efficiencies. The numbers used come from Table V. 

Most systematic uncertainties cancel between / and /, with the exception of charged 
pion and kaon tracking and particle identification. Here, the relevant factor is the charge 
dependence of the efficiencies in data and Monte Carlo simulations. Separate A'^, vr"*", 
and 7r~ tracking and particle ID efficiencies have been determined using the same methods 
that were used to determine overall tracking and particle ID systematic uncertainties. We 
use these efficiencies in data and Monte Carlo to determine systematic errors for the CP 
asymmetries. Kaon tracking produces the largest uncertainty, 0.7% for modes with a charged 
kaon. 

The asymmetries obtained in this analysis are given in Table XI, along with results from 
previous experiments. The uncertainties are of order 1% in all modes, and no mode shows 
evidence of CP violation. Except for the Cabibbo suppressed decay — > K^K~ti^ , our 
results are more precise than previous measurements. We are insensitive to asymmetries 
at the level expected from the Standard Model, the largest of which are a few tenths of a 
percent in modes with a Kg [24] . 
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TABLE XI: The CP asymmetries obtained in this analysis and results from previous experiments. 









CLEO-c 


Previous Results 


Mode 








AcP [/o) 


Reference 




1\ /I 




_n 4. _|_ f) c; _|_ n q 








K-Tr+ 




0.2 ±0.4 ±0.8 


-3.1 ±8.6 


CLEO [25] 




K-TT+ 




0.7 ±0.5 ±0.9 










-7r+ 


-0.5 ±0.4 ±0.9 






D+ ^ 




-tt+ttO 


1.0 ±0.9 ±0.9 






D+ ^ 






-0.6 ± 1.0 ±0.3 


-1.6 ± 1.5 ±0.9 


FOCUS [26] 








0.3 ±0.9 ±0.3 












0.1 ± 1.1 ±0.6 








K+K 


-7r+ 


-0.1 ±1.5 ±0.8 


0.7 ±0.8 


PDG [1] 



XI. CONCLUSIONS 

Using a sample of 281 pb~^ of e+e" DD data obtained with the CLEO-c detector at 
-E'cm = 3.774 GeV, we have measured branching fractions for three hadronic decays and 
six decays. The environment at cc threshold provides a unique opportunity to measure 
these branching fractions. The signals are extremely clean, as illustrated in Figs. 4 and 7, 
and the fact that the double tags are produced without any additional hadrons allows a clean 
determination of the number of produced DD events. In addition, this clean environment 
allows us to directly measure tracking efficiencies, particle identification efficiencies, and 7r° 
reconstruction efficiencies in data. This gives us a good control of systematic uncertainties. 
The single largest systematic uncertainty for the D^ — >■ K~t:^ mode, and several other 
modes, is due to final state radiation. 

The branching fraction results are presented in Table VIII, and the correlation coef- 
ficients among the results are given in Table IX. The branching fractions quoted corre- 
spond to the total inclusive branching fraction including final state radiation photons. Our 
results agree well with (and supersede) our previous measurements based on a 56 pb~^ 
subsample [2] of these data. In all cases the uncertainty of the CLEO-c result reported 
here is less than the uncertainty of the corresponding PDG 2004 [13] average. (We do 
not compare these results to the PDG 2006 [1] averages because the latter include the 
results from the published CLEO-c 56 pb~^ data sample.) Our measurement of the 
reference branching fraction B{D^ K-tt+) = (3.891 ± 0.035 ± 0.059 ± 0.035)% is 
smaller than, but consistent with, that reported recently by the BABAR collaboration [27] , 
BiD'^ K-71+) = (4.007 ±0.037 ±0.070)%. Our result for the reference branching fraction 
bId+ -> K-n+n+) = (9.14 ± 0.10 ± 0.16 ± 0.07)% is substantially more precise than the 
PDG 2004 [13] average. The third errors quoted for our reference branching fractions are 
the systematic errors in our estimates of the effect of final state radiation. Had we not in- 
cluded FSR in our simulations, our quoted branching fractions would have been lower than 
we report; the difference is mode-dependent, ranging from 0.5% to 3% for the branching 
fractions that we measure. 

Our measurements of the production cross sections a{D^D^) — 3.66 ± 0.03 ± 0.06 nb, 
a{D+D-) = 2.91 ±0.03 ±0.05 nb, and a{DD) = 6.57 ±0.04 ±0.10 nb are in good agreement 
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with our earlier measurements using the 56 pb ^ subsample [2] of these data. Again, the 
results reported here supersede the previous measurements. These cross sections agree well 
with the cross sections a{D^W) = 3.39±0.13±0.41 nb and a{D+D-) = 2.68±0.10±0.45 nb 
obtained by combining BES measurements [28] of the branching fractions B{'ijj{3770) — > 
DO^O) = (46.7 ±4.7 ±2.3)% and i3(V'(3770) ^ D+D-) = (36.9 ± 3.7 ± 2.8)%, respectively, 
with the BES measurement [17] of the observed cross section a(e+e~ ■?/'(3770)) = 7.25 ± 
0.27 ± 0.34 nb. Furthermore, our value of the ratio a{e~^e~ — > D^) / a{e^e~ — > D^D^) = 
0.79 ± 0.01 ± 0.01 agrees well with the value a{e+e- D+D-)/a{e+e- D^W) = 
0.79 ± 0.07 ± 0.05 reported by BES [28]. 
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APPENDIX A: SIGNAL Mbc SHAPES 

In this section we describe the form we use for the signal peak in the fits to the Mbc 
distributions for the extraction of signal yields. 

There are four major contributions to the signal line shape. The first is due to the beam 
energy spread. When CESR-c is operating at = 3.774 GeV, the spread in ccntcr-of- 
mass energy is = 2.1 MeV, which is much smaller than the width of the ■?/'(3770). 
The second arises from the effects of initial state radiation (ISR), which reduces slightly the 
center-of-mass energy of the e+e" collision, and — as mentioned in Sec. VI A — produces 
a radiative tail toward larger values of Mbc- The third contribution is the ■0(3770) natural 
line shape, and the fourth contribution is the momentum resolution of the reconstructed D 
candidates. 

The distribution function f^{E) of the '?/'(3770) energy depends on the energy spectra of 
the CESR beams and ISR photons, as well as the ^^^(3770) hue shape. The distribution of 
the total energy E^ot of the e'^e~ pair before ISR is modeled by a single Gaussian, 

9E{Etot; £^cm, <7e) = ^^e-(^--^^-)'/(2'^-), (Al) 

V27r(TE 



where £^cm = '^Eq is the mean total energy of the CESR beams and cte is the energy spread. 
The distribution^ of the energy of ISR photons is taken to be [29] 



where 

/3 = ^ 



h{E,) = E^-\ (A2) 

(A3) 



TT 



E 

2ln' 



rUf. 



* For simplicity, these distribution functions are not normalized to 1. The RooFit [19] fitting package takes 
care of the overall normalization of the distribution functions used in fits. 
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At the i/j{3770) resonance, f3 ^ 0.078. The energy distribution fe+e-{E) of the and e~ 
when they colhde is obtained from an integration of the beam energy spread and the ISR 
photon energy distribution, 

roo 

f,+,-{E)^ h{E^)gE{E + E^)dE^. (A4) 
Jo 

Although the energies of the e+ and beams in CESR are equal, the center-of-mass frame 
of the e^e~ collision and the laboratory frame are slightly different for two reasons. First, 
the beams in CESR-c approach the interaction region at a small crossing angle, 9c ~ 2.5 
mrad, which results in a small boost perpendicular to the axis of the drift chamber system. 
Second, the e+e~ pair is boosted from recoil against ISR photons, whose average momentum 
is of order 3 MeV/c. Monte Carlo simulations demonstrate that the effects of these two 
Lorentz transformations are modest and are readily absorbed in momentum resolution effects 
described below. Hence, we treat the e+e~ center-of-mass frame to be the same as the 
laboratory frame. 

In this analysis, the ■0(3770) natural line shape is taken to be 

^""^^^^ ^ {E^ - MIY + {M^Vt{E)Y' ^^^^ 

where is the mass of the ■?/'(3770). The total width Tt{E) is the sum of the partial 
widths for neutral and charged DD pairs, Vt{E) = ro(£') -\- ^+{E). The numerator V{E) 
is either ro(-E') or V^{E) depending on whether or D^D" events are being fit. These 

partial widths are 



r„(B) = r,e„-3-^^^and 



respectively. In these expressions, is the measured width of the V'(3770), Bo{B^) is 
the branching fraction for the decay of the ■0(3770) to D^D^{D'^D~) pairs, (g+) is the 
momentum of a D^{D~^) of energy E/2, and qoMiQ+Ai) is the momentum of a D^{D^) of 
energy The branching fractions that we used are Bq = 0.57 and B+ = 0.43. The 

parameter r is the Blatt-Weisskopf interaction radius. We use r = 12.3 GeV~^ = 2.4 fm, 
the value favored by our data given the BES mass and width parameters. 

The energy distribution of the ^0(3770) mesons that are produced is obtained by multi- 
plying the e'^e~ energy distribution fe+e-{E) with the cross-section for ^0(3770) production, 

/•oo 

U{E) = fBw{E) / h{E^)gE{E + E^) dE,. (A7) 

J 

The ^0(3770) energy E is related to q, the magnitude of the momentum of the produced 
D and D,hy E ^ 2sjq^c^ + mfjcl Hence, f^{E) can be transformed into a distribution 
function uoiq) for the D momentum. 



uoiq) = U{E) 



dE 



dq 



(A8) 
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The measured D momentum p differs from q due to detector resolution and the effects of 
the two Lorentz transformations relating the center-of-mass frame of the e^e~ collision to the 
laboratory frame. Monte Carlo simulations show that the resulting resolution distribution 
is described well by the sum of three-dimensional Gaussian resolution functions. Each term 
in this sum is given by 

where q is the momentum of the D meson, p is the reconstructed momentum, and is 

the momentum resolution, assumed to be the same for both longitudinal and transverse 
components of p relative to the direction of q. The D is reconstructed from multiple final- 
state particles, and the vector sum of their momenta tends to average out any directional 
dependence. In our fits, there are three such terms, each with a different value of Up (see 
Eq. (7)). In the discussion below, we consider smearing with a single Gaussian, for simplicity. 
The extension to the sum of three Gaussians is straightforward. 

Since the line shape distribution ud^q) depends only on the magnitude g = |q| of 
the D meson momentum, we reduce the three-dimensional momentum resolution function 
g'p(p; q, CTp) to a one-dimensional resolution function r{p;q,ap) for the probability distribu- 
tion of the measured value oip = |p| given the produced value of q. This requires integrating 

9p{Pl CTp) dpdfl over angles transverse to q. In this expression, dpdfl is the the usual 
spherical coordinate volume element and the polar and azimuthal angles of dfl are relative 
to the vector q. Therefore, 

r(p; q, a^) = p" j ^^(p; q, a^) dQ = fe-^^-^)'/^^.^ _ g-(p+,)V(2.|)l _ ^^10) 



The distribution of the reconstructed D momentum, vd{p), is then determined by smearing 
the distribution of the true D momentum, uoiq) of Eq. (A8), with r(p; q, ap), 

roc roo 

Vd{p) = r{p;q,ap)uD{q)dq= r{p;q{E),ap)f^{E) dE. (All) 

Jo J2mD 

Since the measured value of Mbc is a function of the reconstructed momentum p (Eq. (6)), 
the distribution function wd{M-qc) of Mbc is related to voip) by 



wd(Mbc) 



dp 



dMBc 



]\/f roc 

Vd{p) = / r{p- q{E), ap)U{E) dE. (A12) 

P J2mD 



The distribution for double tags, i.e. for Mbc = Mbc(-D) and Mbc = ^bc(-D), is sim- 
ilar to the form developed above for a single Mbc distribution. Since both D mesons are 
produced with the same momentum q, Eq. (All) generalizes to the following probability dis- 
tribution for reconstructing the DD pair with measured momenta p and p given resolutions 
(7p and CTp, 

Vdd{p^p) = / '^(P; f^pyiP; q, ^p)uD{q) dq. (A13) 
Written in terms of Mbc and Mbc, we have 

Mun Mr^ r°° 

«;^^(Mbc, Mbc) = / r(p; q{E), ap)r{p; q{E), ap)U{E) dE. (A14) 

P p J2m£, 

The single and double tag distributions in Mbc cannot be evaluated in a closed form. In 
the fitter, based on the RooFit [19] package, the integrals are implemented numerically. 
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APPENDIX B: SYSTEMATIC UNCERTAINTIES IN CHARGED TRACK, 
AND 7r° RECONSTRUCTION EFFICIENCIES 



We use Monte Carlo simulations to estimate our efficiencies for reconstructing D decays. 
For precision measurements, we must also understand the accuracy with which the Monte 
Carlo events simulate these efficiencies. To determine systematic uncertainties for tt^, K^, 
K^, and 7r° reconstruction efficiencies, we mcasiu^c efficiencies for each particle type in data 
and Monte Carlo simulations using a partial reconstruction technique. We then determine 
the difference, eMc/^data — 1, for tt^, K"^, Kg, and 7r° reconstruction efficiencies, where euc 
is an efficiency found in Monte Carlo simulations and edata is the corresponding efficiency 
found in data. 

We first reconstruct all particles in each event except the particle (X) whose efficiency 
we wish to measure. We calculate the missing mass squared (M^^gg) from the reconstructed 
particles. This variable peaks at M^, the square of the mass of the missing particle X. Then 
we look for the missing particle in each event and separate events into two classes, those for 
which the missing particle was found and those for which it was not found. Peaks in these 
two M^jgg distributions at Mx give the number of times we did and did not find X. From 
these numbers, we calculate the efficiency. This procedure is performed independently with 
data and Monte Carlo samples. 

We use iIj{3770) DD events to measure tt"^ tracking, K"^ tracking, and reconstruc- 
tion efficiencies, and ipi'^S) — > J/ipTiTT events to measure the efficiencies for low- momentum 
TT^s and 7r°s. In DD events, we reconstruct a tag D and all but one of the decay products of 
the D, form M^j^^, and then search for the missing particle. In ip{2S) events, we reconstruct 
J /%p and one of the pions, form M^j^^, and search for the missing 7r+ or tt*'. 

In DD events, we select D candidates using the same selection requirements found in 
Sec. IV, except that we use more restrictive Mbc and A£' requirements: |Mbc — Af^l < 
0.005 GeV/c^ and |A£'| < 0.025 GeV. These requirements produce a clean sample of tags. 
We reconstruct D^ tags in the modes D^ K^-n^ , D'^ — >• ir+7r~7r°, and D^ K'^7r~Ti~^7i~ , 
and we reconstruct D^ tags in the modes D^ K~^7r~7r~ and D~ Kg7r~. In some 
cross-checks we also consider additional tag modes. 

We measure the reconstruction efficiencies for charged pions and kaons in the decays 
D° — s> K~7r~^, D^ — > K~7r'^7r^, and D~^ K~7r~^7r~^, and for mesons in the decay 
D'^ Kg7r~^7r~. In each case, we combine a D^ or D^ with the other particles in the D° or 
D~^ decay. These particles are subject to the selection requirements found in Sec. IV. 

For each of these combinations we calculate the missing mass squared, 

^miss = (Ptot -Pd- Pother)^ (Bl) 

where is the four-momentum of the reconstructed D, Pother is the four-momentum of the 
other particles that were combined with the tag D, and ptot is the four-momentum of the 
e~^e~ pair. In the missing mass squared calculation, we constrain the beam constrained mass 
Mbc of the tag D to the known D mass. That is, we rescale its momentum magnitude to the 
expected value given the beam energy, but leave its direction unchanged. This constraint 
improves the M^jgg resolution. 

For tracking efficiency measurements, we ignore candidates in which the missing momen- 
tum fails a polar angle requirement, | cos^^| < 0.9. This eliminates candidates in which we 
expect the missing particle to be outside of the tracking fiducial volume, where we would not 
be able to detect it. The requirement of | cos^j < 0.9 on the missing momentum is narrower 
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than the angular acceptance of the CLEO-c detector, | cos6\ < 0.93. We choose a narrower 
requirement because the missing momentum direction, determined from the other particles 
in the event, may differ slightly from the true momentum direction of the missing particle. 
We later add an additional uncertainty for tracks that are eliminated by this requirement 
on the missing momentum, but lie within the tracking fiducial volume. 

We next consider all remaining tracks or Kg candidates that pass the requirements found 
in Sec. IV, except we ignore particle identification requirements. If we find a particle that 
forms a good D candidate when combined with the other D decay products, then we have 
found the missing particle. The requirements for a good D candidate are \Mbc — Md\ < 
0.01 GeV/c^ and \A.E\ < 0.05 GeV/c. If we do not find a good D candidate, then we have 
not found the missing particle. 

In events in which the missing particle was found, we fit the clean M^jgg peak at 
with a sum of two Gaussians. A small flat background term is also included in the flts. 

Events in which the particle was not found arc of two types: "inefficient" events in which 
the missing particle was present but not detected, and background events in which it was 
not present at all. The inefficient events form a peak at the particle mass squared with the 
same shape as the peak in events where the particle was found. Therefore, the shapes and 
positions of the inefficient peaks are fixed to match those of the efficient peaks. The fits also 
include one or more terms for the backgrounds. The shapes of the background distributions 
are different in different modes. Fit parameters in data and Monte Carlo events are always 
independent of each other, except for any background shapes that are determined from 
signal Monte Carlo events. 

For TT^ tracking efficiencies, wc make a total of seven measurements of eMc/cdata ~ 1; "the 
efficiency discrepancy between data and Monte Carlo simulations. Three modes are used 
— D'^ —> K~t:~^, —>■ K~7r~^7:^, and —>■ K~7i^n^ — and in the three-body modes we 
measure efficiencies in three momentum bins: 0.2 < p-,^ < 0.5 GeV/c, 0.5 < Ptt < 0.7 GeV/c, 
and > 0.7 GeV/c. For kaons, we make another seven measurements with the same 
modes and momentum bins. We expect that the dependence of efficiency on momentum and 
particle type will be wcll-modclcd in the simulation because efficiency is mostly determined 
by the probability of decay inside the drift chamber, and this probability is expected to be 
well-understood. Therefore, we expect that the discrepancy between data and Monte Carlo 
efficiencies, if any, is similar for pions and kaons and for different momentum ranges. To 
estimate the systematic uncertainty for track reconstruction, we average all measurements 
of eMc/^data — 1 and then add additional uncertainties for the small fraction of tracks not 
included in these measurements (low momentum or high | cos 6*1). 

For Kg reconstruction efficiency, we determine the systematic with the mode D° — > 

1. Charged Pion Tracking Efficiencies 

We illustrate the technique for pion tracking efficiency measurements by describing one 
measurement in detail, and then we present results for all measurements. 

Figure 8 shows missing mass squared plots in which the missing particle is a tt"'" in 
K'TX^-K^ . For data and Monte Carlo events, wc show the distributions when the 
missing tt"'" track was or was not found. In these plots we have combined all three momentum 
bins. 

There are several backgrounds in the plots where no tt"*" track was found. — > 
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FIG. 8: Histograms of and fits to M^j^g distributions from — > K^tt^tt^ decays to determine 
the charged pion efficiency for > 0.2 GeV/c. Figures (a) and (c) are from events in data, and 
(b) and (d) are from events in Monte Carlo simulation. Figures (a) and (b) are from decays in 
which the pion was found, while (c) and (d) are from decays in which the pion was not found. 
The solid curves are fits to the data or Monte Carlo sample; the dashed curves in (c) and (d) are 
background contributions. 

i^'^7r+7r+7r° appears to the right of the signal peak, and the semileptonic decays — »• 
K^Ti^e^Ue and —>■ K^tt^ fi^u^ turn on under the signal peak. All three backgrounds are 
modeled by error functions; this shape accounts for a kinematic threshold with smearing. 
Parameters for these error functions are determined in separate Monte Carlo simulations. 
In the fits to data and generic Monte Carlo events, all yields are allowed to float separately, 
except that we fix the ratio of the backgrounds K~TT^e^Ve and K~tt^ ii^u^j, 

according to their relative efficiencies. The yields and efficiencies in separate momentum 
bins are shown in Table XII. 

The semileptonic backgrounds turn on under the signal peak, so we are sensitive to their 
shape. We determine systematic uncertainties from these backgrounds by varying the widths 
and positions of the two error functions in data. These systematic uncertainties are small 
compared to the statistical uncertainty. The results from this mode are consistent with zero 
efficiency difference between data and Monte Carlo simulations. 

In total, we have seven independent measurements of the differences in charged pion 
tracking efficiencies between data and Monte Carlo simulations — one from — > K~7c^ 
and three each from K^tt^tt^ and K^-k^-k^ . The latter two modes provide 

measurements in each of the three momentum bins. The seven measurements are shown in 
Table XIII. All measurements are consistent with zero difference between data and Monte 
Carlo simulations. The average of these measurements is (0.02 ± 0.26)%, also consistent 
with zero difference. 

We have also measured the pion tracking efficiency from the low-momentum pions in 
ipi'^S) — » J/iIjit~^tt~ with a similar technique. This analysis finds agreement between data 
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TABLE XII: Charged pion yields and tracking efficiencies from — > K~ti^'k^ in three momen- 
tum bins. The systematic uncertainties in the efficiency differences come from varying the shape 
of the — if^7r+^+f£ background in data. 



0.2 < < 0.5 GeV/c 


Data 


Monte Carlo 


Number found 


2766 ± 53 


23174 ± 153 


Number not found 


99 ±13 


875 ± 40 


Efficiency (%) 


96.54 ± 0.44 


96.36 + 0.16 


CMc/edata - 1 (%) 


-0.19 ± 0.49 ± 0.05 (-0.4(7) 


0.5 < < 0.7 GeV/c 


Data 


Monte Carlo 


Number found 


4143 ± 65 


38087 ± 200 


Number not found 


132 ± 16 


990 ± 45 


Efficiency (%) 


96.91 ± 0.37 


97.47 + 0.11 


CRIC/Mala - i (%) 


+0.57 ±0.10 ±0.09 (±1.1(7) 


p^+ > 0.7 GeV/c 


Data 


Monte Carlo 


Number found 


1694 ± 43 


14480 + 125 


Number not found 


47± 14 


345 + 38 


Efficiency (%) 


97.30 ± 0.79 


97.67 + 0.25 


CMc/edata "1 (%) 


+0.38 ± 0.85 ± 0.08 (±0.4(7) 



TABLE XIII: Measurements of the charged pion tracking efficiency differences between data and 
Monte Carlo simulations, and averages of these measurements. In this table, statistical and sys- 
tematic uncertainties are combined. 





CMc/edata - 1 (%) EMc/edata " 1 (%) 


Average 

CMc/edata - 1 (%) 


0.2 < p^+ < 0.5 GeV/c 
0.5 < p^+ < 0.7 GeV/c 
p^+ > 0.7 GeV/c 


-0.32 + 1.34 -0.19 + 0.49 

-1.03 + 2.24 +0.57 + 0.41 
+0.59 + 3.63 +0.38 + 0.85 


-0.21 + 0.46 
+0.52 + 0.40 
+0.39 + 0.83 


L>o K-Tr+ 




-1.25 + 0.71 


Overall average 




+0.02 + 0.26 


and Monte Carlo simulations at the 0.2% level. We do not use this result when computing 
the tracking systematic uncertainties, but it serves to validate that no correction is needed. 

These measurements are combined with the charged kaon tracking efficiency measure- 
ments, described below, to obtain a final tracking efficiency systematic uncertainty. 



2. Charged Kaon Tracking EfRciencies 

We also use the same three D decay modes to measure the charged kaon tracking effi- 
ciency. The procedure is the same as for measuring pion tracking efficiency, except that the 
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missing particle is i^^, and the backgrounds are different. We show the measurement using 
— ^ K~7r^ and then quote results for the other two modes. 
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FIG. 9: Histograms of and fits to Af^j^^ distributions from K^tt^ decays to determine the 

charged kaon efficiency. Figures (a) and (c) are from events in data, and (b) and (d) are from 
events in Monte Carlo simulation. Figures (a) and (b) are from decays in which the kaon was 
found, while (c) and (d) are from decays in which the kaon was not found. The solid curves are fits 
to the data or Monte Carlo sample; the dashed curves in (c) and (d) are background contributions. 

Figure 9 shows plots of M^j^^ for this mode. As in the pion efficiency measurements, 
when the missing kaon is found we see a clean peak. When the kaon track is not found, we 
observe a significant background. This background arises from the decay modes tt+tt^, 
D° K~TT~^Tr^, — ^ TT^/i+z/^, —>■ K~fi~^i'fj_, and other small contributions. The first 
background, tt+tt^, produces a peak at the pion mass squared. Backgrounds from the 

decays — > K^tt^tt^ and —>■ K~ ^'^v^ are modeled as error functions, with parameters 
determined by separate Monte Carlo simulations. We fix the ratio of the —>■ K^-k^t^^ 
and —>■ K^fi'^u^^ yields in our fits, based on the efficiencies and branching fractions of 
these modes. The last background shape is a straight line rising from a cutoff; this shape 
accounts for — > Tr"/!"*"!/^ and other small backgrounds. 

The yields and efficiencies from the fits are shown in Table XIV. The efficiencies in data 
and Monte Carlo simulations agree well. 

In total, we have seven independent measurements of the differences in charged kaon 
tracking efficiencies between data and Monte Carlo simulations — one from — >■ K~tt^ 
and three each from — > K~tt'^tt^ and K~7t~^tt'^. The latter two modes provide 

measurements in each of the three momentum bins. The seven measurements are shown in 
Table XV. All measurements are consistent with zero difference between data and Monte 
Carlo simulations. The average of these measurements is (0.02 ± 0.40)%, also consistent 
with zero difference. 
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TABLE XIV: Charged kaon yields and tracking efficiencies from D — > K tt"*" decays. 



Data Monte Carlo 
Number found 6126 ± 79 59998 ± 248 

Number not found 620 ± 34 5978 ± 99 
Efficiency (%) 90.81 ± 0.47 90.94 ±0.14 

CMc/edata - 1 (%) 0.14 ±0.54 (+0.3a) 



TABLE XV: Measurements of the charged kaon tracking efficiency differences between data and 
Monte Carlo simulations, and averages of these measurements. In this table, statistical and sys- 
tematic uncertainties are combined. 





£,0 ^ ^-^+^0 


D+ K-Tr+Tr+ 


Average 




eMcAdata- 1 (%) 


CMc/edata - 1 (%) 


CMc/edata - 1 (%) 


0.2 < pk- < 0.5 GeV/c 


+1.64 ±2.31 


-2.00+1.20 


-1.23 + 1.06 


0.5 < Pk- < 0.7 GeV/c 


-0.78 ± 1.69 


+1.22 + 1.40 


+0.41 + 1.08 


Pk- > 0.7 GeV/c 


+1.04 ± 1.55 


-0.06+1.26 


+0.38 + 0.98 


K--K+ 






+0.14 + 0.54 


Overall average 






+0.02 + 0.40 



3. Other Considerations and Conclusions for Charged Tracking Efficiencies 

We have measured tracking efficiency systematic uncertainties for both pions and kaons. 
We measure the difference between the Monte Carlo and data efficiencies to be (0.02+0.26)% 
for pions and (0.02 ± 0.40)% for kaons. We expect the pion and kaon tracking efficiencies to 
be highly correlated, and we average the pion and kaon results to obtain an overall average 
for CMc/^data — 1- The average is (0.02 ± 0.22)%. Based on this result, we see no need 
to apply a correction to the Monte Carlo tracking efficiency. We next consider additional 
uncertainties as well as additional cross-checks. 

In the tracking efficiency measurements, we have ignored combinations for which the polar 
angle of the missing track is such that | cos ^| > 0.9. Since tracks reconstructed in the angular 
range 0.90 < | cos^| < 0.93 are, however, used in the branching fraction analysis, differences 
between data and Monte Carlo efficiencies in this region must be considered. We use two 
methods to estimate the magnitude of this effect on the overall tracking efficiency difference 
for the entire angular acceptance (|cos^| < 0.93). First, wc measure tracking efficiency 
in combinations where the missing momentum vector points at an angle |cos^| > 0.9. 
Second, we compare the cos 9 distributions in data and Monte Carlo events for reconstructed 
tracks in D candidates. Both of these methods suggest that the possible effect of tracks 
near the boundary of the active tracking volume on the overall data-Monte Carlo efficiency 
difference is less than 0.2%. Therefore, we add an additional systematic uncertainty of 0.2% 
in quadrature with the other uncertainties on the average difference eMc/^data — 1- 

We have also ignored the lowest-momentum tracks, in particular curlers — tracks whose 
transverse momentum is too low to reach the outer wall of the drift chamber. In each of 
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the D decays whose branching fractions we measure, less than 5% of tracks are curlers. 
We measure tracking efficiency for pion curlers from K^n^n^ and find agreement 

between data and Monte Carlo simulations with a precision better than 2%. A conservative 
upper bound of the effect of curlers on the overall average of eMc/^data — 1 is 5%x2% = 0.1%. 
We add this uncertainty of 0.1% in quadrature with the other uncertainties on eMc/^data — 1- 

We have seen excellent agreement in DD events between data and Monte Carlo tracking 
efficiencies. Similar studies of ip{2S) —>■ J/'^tt+tt" find agreement for both low- momentum 
pions and high-momentum muons with a precision of 0.2%. All of these results indicate 
that no correction to the Monte Carlo tracking efficiency is necessary. To obtain a track- 
ing efficiency systematic uncertainty, we add in quadrature the measured uncertainty on 
CMc/^data — 1, 0.22%, aloug with the additional systematic errors of 0.2% and 0.1%. This 
gives a systematic uncertainty of 0.3% per track, correlated among all tracks. 

We also performed a number of cross-checks to verify consistency between data and Monte 
Carlo simulations in the dependence on polar angle, tag D decay mode, and charge. All 
cross-checks showed good agreement except for the dependence of kaon tracking efficiency 
on kaon charge. We expect that the difference between and K~ tracking efficiencies 
is not larger than a few tenths of a percent, based on hadronic cross sections and the 
amount of material in the beampipe and drift chamber. The efficiencies in Monte Carlo 
simulations agree with this prediction, but the efficiency difference in data between and 
exceeds the difference in Monte Carlo simulations by 1.23 ± 0.61%. Since the average 
efficiency showed no difference between data and Monte Carlo events, this indicates 
that the simulation may have an error of order 0.6%. Therefore, we add an additional 0.6% 
systematic uncertainty for each kaon track, correlated among all kaons. 

4. Kg Reconstruction Efficiencies 

The measurement of Kg — > 7r^7r~ reconstruction efficiencies is similar to the measurement 
of tracking efficiencies. In this case, the goal is to measure the efficiency for correctly 
reconstructing a Kg — > 7r"'"7r~ vertex from a pair of candidate tracks that were found. We 
use D^D^ events in which either D° or decays to Kg'K'^'K~ . We wish to measure the 
efficiency given that the two pions from the K^ decay were found; that is, we wish to 
factor out tracking efficiency from our measurement. Furthermore, we need to eliminate 
— > ir°7r+7r~ and K^g — > 7r°7r° events from our data and Monte Carlo samples. Both of 
these modes would contribute to the peak at Mj^o for not-found candidates, but Kg — > 7r°7r° 
would not contribute to the peak for found candidates and the contribution of K^ ti^tt~ 
to the peak for found candidates would be insignificant. To accomplish both of these goals, 
we require that, in addition to the tag and two pions, the event must contain another 
pair of oppositely charged tracks loosely consistent with the hypothesis of a missing K^g. 
Specifically, the invariant mass must satisfy 0.3 < M(2 tracks) < 0.7 GeV/c^, and the 
magnitude of the vector difference between the pair's momentum and the predicted K^ 
momentum must be less than 60 MeV/c, a value determined from the momentum resolution 
for events in which the K^ was found. Events which do not contain a suitable pair of tracks 
are removed. In the remaining events, we look for a K^ with the standard K^ vertex finder 
using the requirements described in Sec. IV. 

Fake K^ candidates resulting from random combinations of charged pions make the sep- 
aration of M^jgg distributions into signal and background quite comphcated. The overall 
M^igg distribution, before separation into cases where a K^g was or was not found, consists 
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FIG. 10: Histograms of and fits to M^j^g distributions to determine the Kg efficiency. Figures 
(a) and (c) are from events in data, and (b) and (d) are from events in Monte Carlo simulation. 
Figures (a) and (b) are from decays in which the Kg was found, while (c) and (d) are from decays 
in which the Kg was not found. The background peak and deficit are determined by searching 
for Kg candidates in high and low sidebands of the Kg mass. In Figs, (a) and (b), the dashed 
curves are the contributions from fake Kg candidates. In Figs, (c) and (d), the dashed curve is 
the background — a linear function with a deficit due to events in which a fake Kg candidate was 
found — and the solid curve is the total fit function including the signal peak. The area between 
the curves is proportional to the number of Kg mesons not found. 

of a peak from D'^ — > Kg7r~^7T~ events and a non-peaking background from D° ti^ti~t[^tt~ 
events and from —>■ KgTT~^7r~ events in which one or both of the pions from the Kg are 
used in forming M^jgg. However, if the two missing pions in a background event happen to 
have a mass near M^q, corresponding to a missing mass squared near M|-o, they may be 
reconstructed as a fake Kg. As a result, the partitioning of this roughly flat background 
forms a peaking background under the signal peak for events in which a Kg was found, and 
it leaves a corresponding deficit in the background for events in which no Kg candidate was 
found. We estimate the size and shape of this background peak and deficit by searching 
for Kg candidates whose masses lie in high and low sidebands of the Kg mass. We obtain 
separate background estimates for data and Monte Carlo simulations, so our measurements 
are not biased by any discrepancies in the simulation of the background composition. 

Figure 10 shows the M^jgg distributions and fit results. Each signal or background peak 
is fit with the sum of two Gaussians. The background from fake Kg candidates has been 
determined from the Kg mass sidebands. In events where no Kg was found, the background 
is a linear function with a deficit that matches the background peak from fake Kg candidates. 
The deficit in the background is a significant effect, approximately equal in size to the number 
of true Kg mesons that were not found. 

Table XVI shows the yields and the calculated efficiencies. The uncertainties are sta- 
tistical and, where a second uncertainty is listed, systematic. In evaluating the statistical 
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TABLE XVI: Yields and efficiencies for Kg mesons. The statistical uncertainties on the efficient 
and inefficient KgS do not include uncertainty due to the number of fake Kgs; this uncertainty is 
included in evaluating the statistical uncertainty on the efficiency. The systematic uncertainty in 
data comes from widening the background shape. 





Data 


Monte Carlo 


Number of fake KgS 


224 ± 19 


2271 ± 60 


Number of true KgS found 


2754 ± 55 


23, 759 ± 161 


Number of true KgS not found 


143 ± 25 


1564 ± 73 


Efficiency (%) 


95.06 ± 1.06 ± 0.26 


93.82 ± 0.36 


CMc/edata — 1 


-1.30 ± 1.16 ± 0.27% (-1.1(7) 



uncertainty, we have included the uncertainty in the number of fake Kg candidates; this 
affects the numerator of the efficiency but not the denominator. We evaluate systematic un- 
certainty in the shape of the background — specifically, the possibility that the background 
may be wider in data than the Monte Carlo simulation predicts. This systematic uncertainty 
is much smaller than the statistical uncertainty. The efficiencies are high, as expected. In 
fact, much of the inefficiency may be explained by cases where the Kg daughter pions were 
found, but reconstructed poorly. Then they would pass the loose requirement on pairs of 
extra tracks, but not the tighter Kg selection requirements. For example, if one of the pions 
decayed to fiu^, the reconstructed track may have approximately correct momentum, so that 
it passes the loose requirement but fails the Kg vertex finder. 

We obtain our Kg reconstruction systematic uncertainty from the data-Monte Carlo 
difference of —1.30 ± 1.16 ± 0.27%. We have no reason to expect a difference between 
data and Monte Carlo simulations, and the measured discrepancy is consistent with zero. 
Therefore, we make no correction to the Monte Carlo efficiency. We combine the central 
value and uncertainty of the discrepancy in quadrature to obtain a systematic uncertainty of 
1.8%. This systematic error contributes in addition to the tracking systematic uncertainty 
for the two pion tracks. 

5. 7r° Reconstruction Efficiencies 

Using a technique analogous to that used in ip{2S) — ^ J/ipn^ii^ decays, we measure 
the 7r° efficiency in il'{2S) — >• J/ipTT^Ti^ decays. We reconstruct J/ip candidates in the e'^e~ 
and fj.~^fj,~ decay channels. Electron and muon candidates are subject to the charged track 
requirements described in Sec. IV, except that consistency with the pion or kaon hypoth- 
esis is not applied. Electron candidates must also have associated energy deposited in the 
calorimeter approximately equal to the track momentum as well as dE/dx consistent with 
the expectation for electrons. Muons are identified by straw tube chambers embedded in 
the iron return yoke of the superconducting solenoid. Tracks that penetrate to a depth of at 
least three interaction lengths are considered muon candidates. We select J/ip candidates 
from e^e~ and combinations with invariant mass within 50 MeV/c^ of the known J/i/j 

mass [13]. 

We pair these J/ ip candidates with a 7r° candidate satisfying the requirements given in Sec- 
tion IV, and we calculate the M^^^ for the event, which, for ■0(25') — > J/t/jn^n^ decays, peaks 
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FIG. 11: Distributions of ir^ missing mass squared in candidate 'ip{2S) — > J/ip-K^Tr^ events for data 
(points) and Monte Carlo events (histogram). The predicted background level is also shown. The 
vertical arrows demarcate the signal region. Events in which the second vr*^ was found are shown 
in (a) whereas the events where the second vr*^ was not found are shown in (b). 



at M^o- To suppress ip{2S) J/ipn'^ transitions, we also require < 500 MeV/c and p^o < 
500 MeV/c. To further suppress fake 7r° contributions as well as other non-signal ipi^S) de- 
cays (especially ip{2S) J/i>v), we also require (p^o+PmisJ-(P^o-Pmiss)V(0-5 GeVVc^) > 
0.10 GeVVc^ and (p^o + pLs) - (pIo -pLJV(2 GeVVc^) < 0.17 GeVVc^ which selects 
the kinematic region expected to be populated by ip{2S) J/tpn^n^ decays. When the 
event contains a second reconstructed vr^ candidate with M{J/ipn'^n^) — M{J/ip) within 
50 MeV/c^ of the nominal ip{2S)-J /ip mass difference, we consider the 7r° found. 

The resulting M^j^g distributions, separated into the cases where the 7r° was or was not 
found, are shown in Figure 11. For the MC efficiency, we use only events where a J/ipir^Tc^ 
decay is known to have occurred. In data, the non-J/ipir^n'^ contribution is negligible in 
the found-7r° sample and is 2% of the undetected-7r° sample, of which approximately 40% 
comes from J/iprj, 50% from J/il)7r~^7r~ , and 10% from XcJl (primarily Xal) followed by 
Xcj J/'4'l- For both the found- 7r° and the undetected- vr" samples in data, a fake 7r° is 
used to calculate M^iss 6-7% of the entries. These backgrounds peak in the same region 
as the true signal decays, so we obtain the vr" yields by counting the number of entries with 
^miss between —0.05 and 0.10 GeV^/c^ and then subtracting the expected non- J/?/^7r°7r° 
contribution predicted by MC and based on previously measured branching fractions [30]. 
We place a conservative systematic uncertainty of 20% on this subtraction. 

Table XVII gives the overall yields {N{tt^)) and 7r° efficiencies in our MC sample and in 
data. The absolute efficiency difference between data and MC is (—2.60 ± 0.43 ± 0.24)%, 
which corresponds to a relative discrepancy oirj = (edata/cMc) ~ 1 = (—4.37 ±0.72 ±0.41)%. 
The 7r° momentum spectrum in ip{2S) J/ipTT^n^ decays lies below 400 MeV/c, with an 
average momentum of 250 MeV/c. However, in our signal modes -D° — * fC^vr+Tr^, — 
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ii'~7r"'"7r"'"7r°, and K^ti^ti^ , the momentum is typically higher, around 450 MeV/c. 

We extrapolate r] from 250 MeV/c to 450 MeV/c by fitting values of r/ measured in bins of 
Pmiss to a linear function. At 450 MeV/c, rj = (—3.9 ±2.0)%, where the uncertainty includes 
a contribution of 1.8% from the extrapolation. This efficiency correction and systematic 
uncertainty is applied to all vr^s in our analysis. We also examined 77 as a function of 
cos ^miss and found no appreciable dependence on this variable. 

TABLE XVII: Yields and efficiencies for vr'' mesons in data and Monte Carlo samples. 



Sample 


N{tt^) found 


A(7r°) not found 


6,0 (%) 


MC 


86936 ± 295 


59032 ± 243 


59.56 ±0.13 


Data 


8102 ±90 ±3 


6123 ±78 ±61 


56.96 ± 0.41 ± 0.24 



APPENDIX C: INTEGRATED LUMINOSITY DETERMINATION 

In e~^e~ collisions, the most useful final states for measurement of luminosity are e~^e~, 
77, and n'^/J,' because each has a well-known cross section calculable in QED. Each is 
distinctive and not vulnerable to substantial backgrounds. Below we describe the event 
selection criteria and backgrounds as well as the MC simulation used for normalization, 
estimate systematic uncertainties, and finally combine the three normalizations into a single 
integrated luminosity. 

Event selection criteria isolate three classes of events. We require that the number of 
charged particles found in the tracking system with loose track quality requirements must 
be at least two for e~^e~ and /j^/j^ candidates, and must be less than two for 77 candidates. 
For e~^e~ and fi'^fi~ candidate events, each of the two tracks with highest momentum must 
have 0.5 < p/ Eq < 1.1, where p is the momentum of the track and Eq is the beam energy. 
We distinguish muon pair events from Bhabha events using the energy Ec deposited in the 
calorimeter by each of the leptons; this energy is calculated by summing the energies of the 
showers encountered by the helical trajectory of the track. For each muon candidate, we 
require this energy deposit to lie in the range 0.1 < Ec < 0.5 GeV; for Bhabha candidates, 
we require that the ratio E^/p of deposited energy to track momentum must exceed 0.8 for 
one track and 0.5 for the other. The deposited energy requirements for both electrons and 
muons are loose and reject only a small fraction of the signal particles. 

An important signature of these luminosity monitoring events, and therefore a key char- 
acteristic distinguishing them from most potential backgrounds, is that, modulo initial or 
final state radiation, nearly the entire center-of-mass energy should be present in just two 
final state particles. We require that the total energy of the lepton candidate tracks, or two 
highest energy photons for 77 candidates, must exceed 90% of the center-of-mass energy. 
(For leptons, this energy includes recovered bremsstrahlung photons, defined as photon can- 
didates found within 100 mrad of the initial direction of the track; the momenta and energies 
of these photons are then added to the lepton four-momenta.) These requirements accept 
the vast majority of signal events while strongly suppressing backgrounds. 

We require that the two most energetic particles be in the barrel region where material 
in front of the calorimeter is minimized and the detector is hermetic. We require that one 
lepton (or, for 77, photon) make an angle of at least 45° with the beam line, and the other 
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make an angle of at least 40°. (This "one tight, one loose" criterion reduces sensitivity of 
the luminosity to the polar angle measurement or the exact position of the beam collision 
point.) In the 77 final state, the two photons must also be back-to-back in azimuth within 
an acoplanarity angle ^ < 50 mrad. This criterion eliminates essentially all radiative Bhabha 
events with two hard photons that have survived after other 77 criteria have been applied, 
since such events typically have ^ > 150 mrad. 

Cosmic rays in the /x^/i" channel are suppressed by requiring that the tracks be close 
to the measured beam collision point. We calculate the average longitudinal and transverse 
distances of closest approach of the two muon tracks from the collision point and require 
that these distances be less than 4.0 cm in the z direction and 0.1 cm in the x-y plane. We 
determine a small residual background of 1%, estimated to within ~10% of itself, by extrap- 
olating the roughly fiat cosmic ray background from outside to inside the above regions. We 
subtract this background from the \x^\x~ event count. 

We find that trigger efficiencies are essentially 100% by examining events selected with 
independent triggering criteria, i.e.. by using only charged tracks or only calorimeter energy. 

Observed dileptonic cross sections depend not only upon the dominant single photon an- 
nihilation process but also have small contributions from interference with resonance decays. 
For muon pairs at energies near £^cm=3.774 GeV, the effect amounts to -|-0.3% due to the 
J and +0.9% due to the ip{2S). For comparison, the corresponding values are +0.4% and 
—5%, respectively, at Ecm = 3.67 GeV, a continuum point below the ip{2S) where CLEO-c 
has also acquired data. The Bhabha cross section suffers smaller relative effects from this 
interference. We take a systematic uncertainty of 20% of the effect in each case to account 
for possible deviation of these resonances from a pure Breit-Wigner shape so far from their 
peaks. 

There is one other source of non-negligible background for the A*"''/^" final state, namely, 
radiative returns to the '4'{2S) followed by ^(25*) fj.^fi~. In the vicinity of Ecm = 
3.774 GeV, this is estimated to amount to a (0.4 ± 0.1)% background. The analogous 
effect for Bhabha events is only 0.02% due to the large t-channel contribution to the e~^e~ 
cross section; in both cases the estimated background is subtracted from the event count. 

We investigated several other possible backgrounds for any of the three final states and 
found that they contribute at the level of 0.1% or below; these include cross- feed of any of 
the final states into the wrong category of event {e.g., Bhabha events found as 77 or fi~^fi~), 
radiative returns to the J/ip followed hj J/ip ^ radiative returns to ip{2S) followed by 

'0(25') — > XJ/i/j and J/ip — >• i'^i^ , tau-pairs, DD pairs, the direct decay ■0(3770) — >• XJ/ip, 
or a directly-produced 7r"'"7r~ final state. 

We simulated all three final states using the Babayaga [31-33] QED event generator, 
and computed visible cross sections for each after processing events through the detector 
simulation and event selection criteria. At E^ra = 3.774 GeV, the cross sections are 63 nb, 
8.5 nb, and 3.7 nb for the e+e", 77, and /x^/i" final states, respectively. We divide the 
background-subtracted event counts by these cross sections to the determine integrated 
luminosities. Comparisons of polar angle distributions with the respective MC predictions 
are shown in Fig. 12, in which good agreement is observed. The slight excess forward- 
backward asymmetry for data relative to MC in cos is most likely due to interference of 
single photon annihilation with the QED box diagram that has two virtual photons. The box 
diagram is not inchidcd in the Babayaga generator, and neither it nor its interference makes 
a significant contribution to the cross section or measured luminosity when integrated over a 
cos 6*^+ region symmetric about zero. There is also a small systematic difference between the 
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FIG. 12: Distributions of CLEO-c i?cm = 3.774 GeV data (solid circles) and Monte Carlo simula- 
tions (histogram) for the polar angles of the positive lepton (upper two plots) in e^e~ and 
events, respectively, and the mean value of | cos0-y| of the two photons in 77 events. In each case, 
the Monte Carlo histograms are normalized to the numbers of data events. 

data and the Monte Carlo distributions in the mean photon angle (| cos 6.y\) for 77 events. We 
take these small discrepancies into account in determining the Detector Modeling systematic 
errors described in the next paragraph. 

Table XVIII shows the systematic errors assigned for results based on the three final 
states. Detector Modeling errors, including those due to lepton and shower finding and 
reconstruction, dominate, in part due to the natures of electron and photon showers, as 
well as their steep polar angle distributions. We estimate these uncertainties by varying 
selection criteria and from dedicated lepton and photon studies. Integrated luminosity from 
77 {fi~^fi~) events is found to be 2.1% (0.2%) larger than that from Bhabha events; these 
variations are reasonable in light of the systematic errors. Statistical errors are negligible. 
A weighted average of the three values is used for total integrated luminosity, which is 
1.004 times the Bhabha result. Accounting for possible correlations in tracking efficiencies, 
radiative corrections, and interference with direct resonance decays among the three final 
states, we assign a relative uncertainty of 1.0% to the combined integrated luminosity. 
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TABLE XVIII: Summary of systematic errors affecting the luminosity measurements, all in percent. 



Systematic Error (%) 



Source 




77 




Radiative Corrections 


0.5 


1.0 


1.0 


Resonance Interference 


0.1 




0.2 


MC Statistics 


0.1 


0.1 


0.3 


Backgrounds 






0.3 


Trigger 


0.1 


0.1 


0.1 


Detector Modeling 


1.0 


1.0 


0.6 


Total in Quadrature 


1.1 


1.4 


1.3 



In summary, we utilize three QED reactions to measure CLEO-c integrated luminosities, 
and we find that the results are consistent with one another. For the data sample used in 
this analysis of hadronic D decays we find / Cdt — 281.5 ± 2.8 pb~^. 
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